Analysis of singular 

solutions for two .nonlinear 
wave equations 



Candidate Supervisor 
Massimo Fonte Prof. Alberto Bressan 



Thesis submitted for the degree of Doctor Philosophiae 
Academic Year 2004/2005 



Contents 



Introduction 3 

I The Camassa-Holm equation 9 

1 The Camassa-Holm equation 11 

1.1 Non-local formulation 11 

1.1.1 The main results 14 

1.2 The Huntcr-Saxton equation 16 

1.2.1 A transportation map 18 

1.2.2 The stability issue of a system of ODE 19 

2 Conservative multi-peakon solution 21 

2.1 Multi-peakon solution in the real line 22 

2.2 Approximation of the initial data 25 

2.3 Periodic multi-peakon 26 

3 Distance defined by optimal transportation problem 33 

3.1 A distance functional in the spatially periodic case 33 

3.1.1 Continuity in time of the distance functional 39 

3.1.2 Continuity w.r.t. the initial data 44 

3.2 A priori bounds 49 

3.3 Definition of the distance in the real line 52 

3.3.1 Comparison with other topologies 54 

3.3.2 Continuity of solutions w.r.t initial data 56 

3.4 Proof of the main theorems 60 

3.5 Uniqueness 63 

II The discrete Boltzmann equation 67 

4 Symmetry groups of differential equations 69 

4.1 Group and diifcrcntial equations 69 

4.2 Symmetries and blow-up 75 

4.3 A group of symmetry for the discrete Boltzmann equation .... 78 



2 CONTENTS 



5 The two dimensional Broadwell model 81 

5.1 The discrete Boltzmann equation 81 

5.2 Coordinate rescaling 82 

5.3 The two-dimensional Broadwell model 83 

5.4 Refined blow-up estimates 86 

5.5 A tentative blow-up scenario 91 

Bibliography 94 



Introduction 



This thesis deals with two strongly nonlinear evolution Partial Differential Equa- 
tion (in the following named P.D.E.) arising from mathematical physics. The 
first one was introduced first by Fokas and Fuchssteiner [28] as a bi-Hamiltonian 
equation, and then was rediscovered by R. Camassa and D.D. Holm [13] as an 

higher order level of approximation of the unidirectional shallow water wave 
equation than the Korteweg-de Vries equation [38]. It can be written as 

u{t,x) : R X R h-* R 

(1) 

Ut + 2ku 

X ^xxt ~t~ ^UUx — 2UxUxx ~t~ ^'^xxxi 

here the unknown u{t, x) represents the water's free surface over a flat bed and 

K is a constant related to the critical shallow-water wave speed (see also [37] 
for an alternative derivation as an hyperelastic-rod wave equation). We refer 
to this equation as to the Camassa-Holm equation, in honour to the first two 
authors which found a physical meaning stemming from the Euler equation. 

The second PDE we want to study is a system of hyperbolic equations with 
quadratic source 

' (ui,...,MAr)(t,a;) :R+ xR2 h^M^ 

V- (2) 
{ui)t + Ci ■ VxUi = aijkUjUk, for alH = 1 . . . TV 

jk 

which is a discretization of the velocities in the plane R^ for the Boltzmann 
equation 

f{t, x,^) :R+ xR^ xR^ >-^R 

dtfit, a;, + ^ • V,/(t, X, = Q(/, f){t, x, . 

The nonlinear nature of these equations leads to the possibility of blow up 

in finite time either for the solution itself, or for the gradient of the solution. 
The typical situation of blow up in finite time is given by the following Ordinary 
Differential Equation (O.D.E.) 

^v = -v'^, v{0)=vo 
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where V should be cither the sohition u or the gradient u^, possibly computed 
along the characteristic curves of the equation. It is well known that the solution 
of this equation has the behaviour where T depends on the initial data 

and it is the time of blow-up whenever vq < 0. 

This situation can occur for the solutions of the Camassa-Holm equation in 
the limit case k = 0. With this condition, the equation (1) may be rewritten in 
nonlocal form as 

1 

Since the i/^norm is conserved for regular solutions, the L°° norm of u is 
bounded, namely < < Ve. Arguing as in the Steepening Lemma 

(see [14]), in [20, theorem 4.1] the authors prove that smooth solution to (3) may 
not be globally defined. Let start from an odd initial data u G H^{]&) which 
has an inflection point in 0, u{0) < and consider the evolution of the slope at 
the inflection point s{t) = Ux{t, 0) < 0. The computation of the function s gives 
the diff'erential inequality 

If Ux(0) is sufflciently small, we get 

1 1 t 
s(t) - s(0) 2 

and then the slope becomes vertical at flnite time. However, the singularities 

thus can occur only in form of wave breaking (sec also [39, 18]), in fact even if 
its slope can become unbounded at finite time, the solution remains bounded, 
because of the inequality ||u(f)||i:,oo < y/E. 

The aim of Part I of this thesis (see also [10, 29]) is to construct a continuous 
semigroup of global solutions in two main cases: 

1. on the space -ffp^r of spatially periodic functions, locally in H^{R); 

2. on a domain of functions with a certain exponential decay at x — > ±oo. 

Result of existence of solutions can be found in [48, 49], and [17, 16] where the 
authors added a small diffusion term to the right hand side of 3 and obtained 
solution of the original equation as a vanishing viscosity limit. On the other 
hand, in [7] was developed an alternative technique, which relies on a new set 
of dependent and independent variables with the specific purpose to resolve all 
the singularities. With this change of variable the solution can be obtained as 
the unique fixed point of a contractive transformation. In Chapter 2 we present 
yet another approach based on the Hamiltonian structure of the Camassa-Holm 
equation. We shall construct the semigroup of global solution starting from 



-\x\ 
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explicit solutions of the Camassa-Holm equation with initial condition in form 
of multipeakon function 

N 

The motivation of this choice is given by the form of traveling wave solution 
(see [13, 20, 21, 22]). Looking for solution of the equation (3) in the traveling 
wave form u{t,x) = U{x — ct), with a function U that vanishes at infinity, one 
obtains the function U = ce"'^"^*', which is a peaked soliton (from this fact 
derives the shortened term peakon). The multipeakon functions are stable, in 
fact not only a single peakon subject to (3) evolves with this form, but also the 
evolution of a superposition of traveling wave (e.g. initial data like uq) remains 
of the same shape 

N 

«(i,x)=^p,(t)e-l--«^(*)l. 

i=i 

The reader can see also [2] for a recursive reconstruction of the multipeakon 
solutions, and [24] which prove that multipeakon solutions are orbitally stable, 
i.e. stable under a general nature of perturbations. 

In [33] the authors prove the existence of a global multipeakon solution 
when the strengths pi are positive for all ?' = 1 . . . iV. In this case the crucial 
fact is that no interaction between the peakons occurs, and then the gradi- 
ent remains bounded, which yields existence and uniqueness of the coefficients 
Pi{t), . . . ,pN{t) and qi{t),. . . , qwit) that are solutions of the Hamiltonian sys- 
tem 

( ^ 

(4) 

Pi = Pi^Pj sign(gi - q'jOe"!'*"'^!, 

with Hamiltonian H = ^ piPje~l^*~^^ l. 

However, a general initial data contains both positive and negative peakons, 
as in the example of the peakon-antipeakon interaction: one positive peakon 
with strength p, centered in — g, moves forward and one negative anti-peakon in 
q, with strength —p moves backward. The evolution of the system produces the 
overlapping of the two peakons at finite time i = r, so that q (see Figure 
!)• 

There are infinitely many ways to extend the solution after the time t of the 
interaction, for example the vanishing viscosity approach in [17, 16] singles out 
the dissipative solutions. As far as the example of peakon-antipeakon interaction 
is concerned, the vanishing viscosity approach selects the sohition that, after 
the time r, is = 0: all the energy E is lost. In section 2.1 we shall construct 
a conservative solution, i.e. a solution for which the quantity E is constant 
for a.e. time t. At the interaction time the energy E = is described 
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Figure 1: Peakon-antipeakon interaction. 



by a Dirac measure entirely concentrated at the single interaction point. After 
the interaction, a positive and a negative pcakon emerge, whose strengths are 
uniquely determined by imposing the conservation of the total energy. 

In Chapter 3 we shall discuss the issue of the uniqueness and the stability. 

Stability. The multipeakon solutions form a continuous semigroup whose do- 
main is dense cither in (K.) or in H^^^. The main novel feature in our approach 
is the construction of a metric J(-, •) on the space (or -f^pgr) determined by 
an optimal transportation problem. While the semigroup generated by (3) is 
not eve continuous w.r.t. the distance, we show that it is Lipschitz contin- 
uous w.r.t. our new distance functional J. The reader can see [47] for earlier 
applications of distances defined in term of optimal transportation problems, 
and [8] , in which the authors recover a semigroup of dissipative sohition for the 
Hunter-Saxton equation [34] (see also [11] for a fixed-point approach for both 
conservative and dissipative solutions). 

The main wcU-posedness result is provided by a Gronwall-type lemma, stem- 
ming from the inequality 



whenever u and v are two multipeakon solutions (see Section 3.1.2). 

Uniqueness. Example 3.1 in Section 3.5 shows that a solution of 3 need not 
be unique. Roughly speaking, every shifted antipeakon-peakon couple is also 
a conservative solution. A conservative solution can be characterized by an 
additional linear transport equation, accounting for the conservation of the total 
energy. It can be done by the following heuristic idea. 

We can think that the absolutely continuous measure /xt, which satisfies 
d/if = {u^ + u^j.)dC tends, as f ^ t, to a Dirac measure with support in 0. 
We introduce thus a further equation for the measure /Uj, whose absolutely 
continuous part is v? + u^, in the following way. Since whenever + is 
regular it satisfies the equation 



^J{u,v)<C{t)-J{u,v) 




) 



= I{u) 



it suggests that fit provides a measure- valued solution of 

dtlJL + {uiJ.)a: = f{u) . 



(5) 



7 




Figure 2: Two particle interaction 



Our result (Theorem 1.3) shows that every solution (u,/it) of (3)-(5), such that 
the absolutely continuous part of /Uj has density + u^, must coincide with the 
one provided by multipeakon approach. 

Part II of the thesis is devoted to the analysis of blow-up for the discrete 
Boltzmann equation (2). Such a equation is obtained by considering a rarefied 
gas for which is supposed that the particles can move only along a finite number 
of direction characterized by the vectors Ci, . . . ,cjv. The unknowns Ui repre- 
sent the densities of particles which travel at speed Cj. By a collision, a pair 
of incoming particles with speeds 0^, Cj is replaced by a new pair of particles 
say c/c, Ci. The rate at which such collision occur is given by UijkUjUk The con- 
centration Ui is thus increasing (or at least is constant) when interact particles 
of speed different to Cj, decreasing when an i— particle collides with someone 
other. Then, the coefficients a^fc are non negative when j, k ^ i and negative 
when either j = i ot k = i. 

If the initial data is suitably small, the solution remains uniformly bounded 
for all times [4]. For large initial data, on the other hand, the global existence 
and stability of solutions are known only in the one-dimensional case [3, 32, 45]. 
Since the right hand side has quadratic growth, it might happen that the solution 
blows up in finite time. Examples where the L°° norm of the solution becomes 
arbitrarily large as f ^ oo are easy to construct [35] . 

In Chapter 5 (sec [9]) we focus our analysis on the two-dimensional Broadwell 
model (see, for example, [12, 46, 15] for a description of the model) and examine 
the possibility that blow-up actually occurs in finite time. In this model the 
permitted direction are 

Cl = (1, 1), C2 = (1, -1), C3 = (-1, -1), C4 = (-1, 1) 

and the particles have a diamond-shape (see Figure 2). 

As wc will show with the theory developed in Chapter 4, since the equations 
(2) admit a natural symmetry group (sec Section 4.3, and [42] for a more general 
theory) , one can perform an asymptotic rescaling of variables and ask whether 
there is a blow-up solution which, in the rescaled variables, converges to a steady 
state. This technique has been widely used to study blow-up singularities of 
reaction-diffusion equations with superlinear forcing terms [30, 31]. See also 
[36] for an example of self-similar blow-up for hyperbolic conservation laws. 
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Our main results show is an a-priori bound on the blow up rate in the L' 
norm. Namely, if blow-up occurs at time T, then one has 



This means that the blow-up rate must be different from the natural growth 



In the final section of Chapter 5 we discuss a possible scenario for blow- 
up. The analysis highlights how carefully chosen should be the initial data, if 
blow-up is ever to happen. This suggests that finite time blow-up is a highly non- 
generic phenomenon, something one would not expect to encounter in numerical 
simulations. 



Mi)||L=o > 5 




T-t 




which would be obtained in case of a quadratic 
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The Camassa-Holm 
equation 



Chapter 1 



The Camassa-Holm 
equation 



The Camassa-Holm equation 



ut + 2ku. 



'X 



Uxxt + '■iuU. 



'X 



2UxUxx + uu. 



XXX 



arises from a higher order level of approximation of the asymptotic expansion of 
the Euler's equations for a shallow water wave theory. Here we do not enter into 
deep details of the interpretation of such a equation, for the physical motivations 
we refer to [13], [22], [23], [37]. 

In the following we focus our attention and we refer to Camassa-Holm equa- 
tion the previous equation with k; = 0. 

1.1 Non-local formulation 

The Camassa-Holm equation can be written as a scalar conservation law with 
an additional integro- differential term: 



Earlier results on the existence and uniqueness of solutions can be found in [48] , 
[49]. One can regard (1.1) as an evolution equation on a space of absolutely con- 
tinuous functions with derivatives Ux € L"^. In the smooth case, differentiating 
(1.1) w.r.t. X one obtains 



wt + (wV2)x + ^x = 0, 



(1.1) 



where P is defined as a convolution: 




(1.2) 




(1.3) 
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Multiplying (1.1) by u and (1.3) by we obtain the two conservation laws with 
source term 



(1.4) 
(1.5) 



As a consequence, for regular solutions the total energy 
E{t) = j [u'^{t,x) +ul{t,x)\ dx 

remains constant in time. 

As in the case of scalar conservation laws, by the strong nonlinearity of 
the equations, solutions with smooth initial data can lose regularity in finite 
time. For the Camassa-Holm equation (1.1), however, the uniform bound on 
||tta;||i2 guarantees that only the L°° norm of the gradient can blow up, while 
the solution u itself remains Holder continuous at all times. 

In order to construct global in time solutions, two main approaches have 
recently been introduced. On one hand, one can add a small diffusion term in 
the right hand side of (1.1), and recover solutions of the original equations as 
a vanishing viscosity limit [17, 16]. An alternative technique, developed in [7], 
relies on a new set of independent and dependent variables, specifically designed 
with the aim of "resolving" all singularities. In terms of these new variables, 
the solution to the Cauchy problem becomes regular for all times, and can be 
obtained as the unique fixed point of a contractive transformation. 

In the present chapter, we implement yet another approach to the Camassa- 
Holm equation. As a starting point we consider all multi-peakon solutions, of 
the form 



JV 



u{t,x) = '^pi{t)e~ 



\x-gi{t)\ 



(1.6) 



These are obtained by solving the system of O.D.E's 



Pi = '^PiPjsign{qi- qj)e 
It is well known that this can be written in hamiltonian form: 



(1.7) 



d 

Qi = -^H(p, q) , 

d 

Pi = -^H{p,q), 
dqi 



H{p,q) = -J2PiPj^~^'"~'"^ 



1,3 



If all the coefficients pi are initially positive, then they remain positive and 
bounded for all times. The solution u = u{t,x) is thus uniformly Lipschitz 



1.1 Non-local formulation 
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continuous. Wc stress, however, that here we are not making any assumption 
about the signs of the pi. In a typical situation, two peakons can cross at a 
finite time r. As f — > r— their strengths Pi,Pj and positions qi, qj will satisfy 

Pi{t)^ +00, pj{t) ^-oo, Pi{t) +Pj{t) ^ p, (1.8) 

<li{t)-^q, qjit)^q, qi{t) < qj{t) for t < t, (1.9) 

for some p,q G M. Moreover, ||Ma;(i)||^oo ^ oo. In this case, we will show 
that there exists a unique way to extend the multi-peakon solution beyond the 
interaction time, so that the total energy is conserved. 

Having constructed a set of "multi-peakon solutions", our main goal is to 
show that these solutions form a continuous semigroup, whose domain is dense 
in the space ff^(R). Taking the unique continuous extension, we thus obtain a 
continuous semigroup of solutions of (1.1), defined on the entire space H^. 

One easily checks that the flow map : u{0) u{t) cannot be continuous 
as a map from into itself, or from into itself. Distances defined in terms of 
convex norms perform well in connection with linear problems, but occasionally 
fail when nonlinear features become dominant. In the present setting, we con- 
struct a new distance J{u, v) between functions u,v € H^, defined by a problem 
of optimal transportation. Roughly speaking, J{u, v) will be the minimum cost 
in order to transport the mass distribution with density 1 + located on the 
graph of u onto the mass distribution with density l + v"^ located on the graph of 
V. See Section 3.1 for details. With this definition of distance, our main result 
shows that 

p{u{t), v{t)) < C ■ J{u{t), v{t)) 

for some constant C and any couple of multi-peakon solutions u,v. Moreover, 
J{un,u) implies the uniform convergence — 0. The distance 

functional J thus provides the ideal tool to measure continuous dependence on 
the initial data for solutions to the Camassa-Holm equation. Earlier applications 
of distances defined in terms of optimal transportation problems can be found in 
the monograph [47]. The issue of uniqueness of solutions must here be discussed 
in greater detail. For a multi-peakon solution, as long as all coefficients pi remain 
bounded, the solution to the system of ODE's (1.7) is clearly unique. For each 
time t, call Ht the measure having density u'^{t)+u'^{t) w.r.t. Lebesgue measure. 
Consider a time r where a positive and a negative pcakon collide, according to 
(1.8)-(1.9). Ast ^ r— , we have the weak convergence /U* ^ for some positive 
measure which typically contains a Dirac mass at the point q. By energy 
conservation, we thus have 

J [u'^{t,x)+u1{t,x)] dx+Hrdq}) = lira J [u'^{t,x)+u'^{t,x)] dx = E{t—) . 

There are now two natural ways to prolong the multi-peakon solution beyond 
time r: a conservative solution, such that 



E{t) = j [u^{t, x) + ul{t, x)] dx = E{t-) 



t> T , 
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or a dissipative solution, where all the energy concentrated at the point q is lost. 
In this case 



E{t) = j [u'^it, x) + ulit, x)] dx = E{t-) - Mr ({«}) 



t >T. 



For t > T, the dissipative solution is obtained by simply replacing the two 
pcakons Pi,Pj with one single pcakon of strength p, located at x = q. On the 
other hand, as we will show in Section 2.1, the conservative solution contains two 
peakons emerging from the point q. As t — > t+, their strengths and positions 
satisfy again (1.8), while (1.9) is replaced by 

qi{t)^q, qj{t)^q, qi{t) > qj{t) fov t > r . (1.10) 

The vanishing viscosity approach in [17, 16] singles out the dissipative solutions. 
These can also be characterized by the Oleinik type estimate 

u^{t,x) < C{l+t-^), 

valid for f > at a.e. x e K. On the other hand, the coordinate transformation 

approach in [7] and the present one, based on optimal transport metrics, appear 
to be well suited for the study of both conservative and dissipative solutions. 

In the following chapters we focus on conservative solutions to the Camassa- 
Holm equation. Wc start with the study of the spatially periodicity because 
it allows us to concentrate on the heart of the matter, i.e. the uniqueness and 
stability of solutions beyond the time of singularity formation. It will spare us 
some technicalities, such as the analysis of the tail decay of u, £is x ^ ±oo. 
In this respect we shall discuss the decay analysis of solutions in Section 3.2 of 
Chapter 3. 

The main ingredients can already be found in the paper [8], devoted to 
dissipative solutions of the Hunter-Saxton equation. 
As initial data, we take 

u{0,x)=u{x), (1.11) 

with u in the space of absolutely continuous functions u with derivative 
Ux G L^- To fix the ideas, we assume that the period of a spatially periodic 
function in the space -ffpgr is 1, so that 

u{x + l) = u{x) xeR. 

On iJpgj we shall use the norm 

/ /.I 2 /.I 2 \ 

WMIhI^^ = \u{x)\ dx + J \uxix)\ dxj 

1.1.1 The main results 



In this section we state the main results of Part I of this thesis. We shall write 
them for the spatially periodic case. 



1.1 Non-local formulation 
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Theorem 1.1. For each initial data u G -ffpor, there exists a solution u{-) of the 
Cauchy problem (1.1), (1.11). Namely, the mapt i-^ u{t) is Lipschitz continuous 
from M into L'^^j., satisfies (1-11) at time t = Q, and the identity 

= -UU:^ - Px (1.12) 

is satisfied as an equality between elements in L^^^. at a.e. time t £ M.. This 
same map t^ u{t) is continuously differentiable from R into L^gj. and satisfies 
(1.12) at a.e. time t e K, for all p € [1,2[. The above solution is conservative 
in the sense that, for a.e. t G R, 

E{t) = [ [u^{t, x) + ul{t, x)] dx = E^ = [ [tl^ix) + ul{x)] dx . (1.13) 
Jo Jo 

Theorem 1.2. Conservative solutions to (1.1) can be constructed so that they 
constitute a continuous flow Namely, there exists a distance functional J on 
Hpgj. such that 

^■\\u- ^;|Uj„ < J{u, v)<C-\\u- v\\hi^^ (1.14) 

for all u,v G -ffpgj and some constant C uniformly valid on bounded sets of 
iJpci- . Moreover, for any two solutions u{t) = ^tu, v{t) = of (1.1), the 
map 1 1-^ J{u{t), v{t)) satisfies 

J{u{t), u) < Ci • |i| , (1.15) 

J{u{t),v{t)) < J(zi,t;)-e^^l*l (1.16) 

for a.e. t € M. and constants C\,C2, uniformly valid as u,v range on bounded 
sets of Hpg^ . 



The previous results can be extended to the following space of functions 
which exponential decay: let a e]0, 1[, then 

Xc = {ue H'<m.) s.t. / [u'^ix) + u^(a;)]e"l^l 
I Jr 

It is not so restrictive one can think, in fact the pcakon functions (sec Section 
2.1), the natural solitary waves of the Camassa-Holm equation which have the 
soliton properties, belong to it. 

Somewhat surprisingly, all the properties stated in Theorem 1.1 are still 
not strong enough to single out a unique solution. To achieve uniqueness, an 
additional condition is needed. 



■ 
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Theorem 1.3. Conservative solutions 1 1— > u{t) of (1.1) can he constructed with 

the following additional property: 

For each t G R, call fit the absolutely continuous measure having density + 
w.r.t. Lebesgue measure. Then, by possibly redefining /i^ on a set of times of 
measure zero, the map t^ fit is continuous w.r.t. the topology of weak conver- 
gence of measures. It provides a measure-valued solution to the conservation 
law 

Wt + {uw)x = {u^ - 2uP")a. . 

The solution of the Cauchy problem (1.1), (l.H) satisfying the properties 
stated in Theorem 1.1 and this additional condition is unique. 

In Section 2.1 we derive some elementary properties of multi-peakon solu- 
tions and show that any initial data can be approximated in H^^^ by a finite 
sum of peakons. In Section 3.1 we introduce our distance functional J (a, v) and 
study its relations with other distances defined by Sobolev norms. The conti- 
nuity of the flow (1.1), together with the key estimates (1.15)-(1.16) are then 
proved in the following two sections. The proofs of Theorems 1.1 and 1.2 are 
completed in Section 3.4. The uniqueness result stated in Theorem 1.3 is proved 
in Section 3.5. As a corollary, we also show that in a multi-peakon solution the 
only possible interactions involve exactly two peakons: one positive and one 
negative. In particular, no triple interactions can ever occur. 

Now we dedicate the rest of this chapter to exhibit an example which will be 
the start point of the tc;clmique we shall develop in Chapter 3. In particular, we 
shall see how the optimal transportation theory can be useful nonlinear equation 
like Camassa-Holm equation. The key point is to define a Monge-Kantorovich 
like metric for a space of Radon measures. In Section 1.2 wc start from the 
Hunter-Saxton equation [8] for give a brief heuristic idea for how this technique 
is involved in. 

1.2 The Hunter-Saxton equation 

The Hunter-Saxton equation describes the propagation of waves in a massive 
vector field of a nematic liquid crystal. Since the physical interpretation and its 

derivation arc beyond to the description of this thesis, wc refer to [34, 8]. It can 
be written in a non-local formulation as a conservation law with a source term: 

+ (t) ^ i (/I " ID ^ ^^'^^^ 

where 

• f > is the time variable, 



1.2 The Hunter-Saxton equation 
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• a; e M is the space variable in a reference frame, 

• u{t,x) E M. is related to the orientation of the liquid crystal molecules in 

the position x at time t. 

Suppose that there exists a smooth solution to (1.17). To the Hunter-Saxton 
equation we can associate the following two conservation laws 

{ul)t + iuu% = 0, (1.19) 

which are obtained by computing the derivative of the equation (1.17) w.r.t. 
the X variable and then multiplying it by to achieve the second one. 

A further conservation law is satisfied by the source term Q^{t, x). Since for 
smooth solutions (1.19) yields the conservation of the energy 

Eo = E{t) = I ul{t,x)dx 
Js. 

the function can be expressed in the following way 

Qit,x) = -^ + ^J ul{t,y)dy. 
By deriving w.r.t. t we obtain 

Qt + uQ^ = 0. (1.20) 
The function Q" is constant along the characteristic curves ^u{t,y) defined by 

^C«(i, y) = u{t, ut, y)), Uo, y) = y- (i-2i) 

Let us remark that the previous equations holds for all the time t in which u 
is a classical solution. It can be seen by the method of characteristics that if 
uo ^ is a smooth initial data and for some xq we have uox{xo) < 0, along its 
outgoing characteristic the gradient blows up in finite time. Since the quantity 
E{t) remains bounded also at the time of blow up, we can think that a finite 
amount of energy will be concentrated at the point of blow-up. In [8] the authors 
focus their attention on solutions that dissipate this quantity of energy. As far 
as the conservative solution is concerned, equation (1.19) will be satisfied in 
sense of measures, i.e. it means that thinking at a measure nt with absolutely 
continuous part which satisfies d/it = u'^{t,-)djC (here with C we indicate the 
Lebesgue measure), it satisfies 

dtUt + da;{ufit) = inV'. 

To find a conservative solution to the Cauchy problem (1.17) with finite energy 
smooth initial condition uq means then to find a couple {u{t), ^t) which satisfies 



18 



The Camassa-Holm equation 



the following system of conservation laws: let /ig be the absolutely continuous 
measure w.r.t. Lcbesgue measure defined by dfxo = uo^dC, then 

/u^\ 1 1 
9tu + 5^ (yj = --Mt(M) + -Mt(] - oo, a;]) u{0,x) = u{x), ^2) 

dtl^t + dx{ufit) = Mtlt=o = Mo- 

Due to the nonlinearity of the problem, as shows [8, Example 2] for the dis- 
sipative solution of the Hunter-Saxton equation, we can aspect that the usual 
"strong" distance stemming from convex norm is not useful in order to construct 
a continuous semigroup of solutions. In the following section, we give a sketch 
of the construction of a metric which yields continuity of solution with respect 
to the initial data. 



1.2.1 A transportation map 

Let uq and vq be two initial data whose associated measures /iQ and /iQ have 
the same total mass AtJ(R) = A*o(^)- Suppose that such a initial data are not 
constant in any interval of M, so that the functions and Q'"° are absolutely 
continuous and increasing. 




Figure 1.1: The map 



We can thus define a continuous map ^0 for which at every a; € R it asso- 
ciates the unique point \E'o(a:) such that 

Q^°{x)=Q^°iMx)) (1.23) 

(see Figure 1.1). Let us Remark that the \l/o is an increasing function, in fact 
since and Q'"° are increasing, then x < y implies 

Q^^i^oix)) = g"°(x) < Q"«(2/) = Q"°(*o(2/)), 

so *o(a;) < *o(l/) holds. 

Now we want to see how the map evolves in time. Since by (1.20) Q", 
are conserved along the caracteristic curves (1-21), the equalities 

Q"(t,^„(t,a;)) =Q^''(x), Q"(i,C„(i,x)) = g''°(x) for all a; e R 
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yields the definition of the transportation map 

M^u{t,X))=^y{t, M^))- 



(1.24) 



1.2.2 The stability issue of a system of ODE 

As in [7, Section 3], we compute a change of variables in order to obtain a system 
of ODE with Lipschitz vector field. Let suppose that the initial data uq is in 
H'^{R). Let set 

oj = 2 arctan(w2;), 

then CO belong to the unit circle T = [0, 2tt], with and 27r identified. Computing 
the derivative of u along the characteristics, having in mind the equation (1.18) 
for Ux we have 



1 



-tan2(w/2) 
1 + tan2(w/2) 



It is thus natural to consider the new unknowns which take values in the space 
R2 X T, 

x" = x"(i,y)= uit,Ut,y)) 
V '^u(i,C«(*,y)) 

and write the corresponding Cauchy problem 



QHt,Ut,y)) I =/(x«(t,y)) 
uJu{t,^u{t,y)) 



with initial data 



X"(0,y) 



(1.25) 



(1.26) 



uo{y) 

2arctan(uox(y)) 

We remark that by definition, Q" is far from to be Lipschitz continuous, then 
we cannot suppose that in the previous system of ODE the function / is a 
Lipschitz vector field. Hence, to overcome this lack of Lipschitz continuity, we 
shall make use of the function '^t introduced in 1.2.1. Let X" and X'^ be two 
solution of the Cauchy problem (1.25)-(1.26), corresponding to the initial data 
uo, vo respectively. We allow ourselves to make an abuse of notation by defining 
the map in the following way: 

*t(X"(t,y))=X''(t,*o(y)). 

We gain a sort of Lipschitz continuity for the function / if we restrict it on 
the manifold located by Let us compute the difference of the vector field / 
evaluated in the points X.'^{t,y) and '3>j(X"(t, y)). 



|/(X"(i,j/))-/(*,(X"(i,2/)))| 



u{t,£,u{t,y)) 

Q''{t,Ut,y))' 
i^u{t,^u{t,y)) 



■v{t,Ut,'^o{y))) 

Q^{t,ut,My))) 
uy{t,^y{t, ^o{y))) 
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since by definition of the map and by the equation (1.20) we have the identity 

Qn*, Ut, y)) - Q'it, Ut, *o(y))) = o 

we deduce the following estimate for the vector field / 

|/(X"(t,2/))-/(vl/t(X«(f,2/)))| < |X"(t,y)-X"(t,y)|. 

From the previous inequality we can prove that the difference of the two solutions 
X"{t, y) and X'"{t, "^oiv)) can be estimated by the initial data. In fact, Gronwall 
Lemma applied to the inequality 

\X.^{t,y)-X.^{t,My))\ < |X«(0,y)-X''(0,*o(2/))| 

+ / |/(X«(s,y))-/(X"(s,vl/o(y)))|ds 
Jo 

yields the estimate 

|X"(t,y) - x-{t,^o{y))\ < e^lX^oCj/) - x^'^{My))\- 

The previous inequality suggests how to introduce a new distance in order to 
obtain a stability result for solutions of the Hunter-Saxton equation. For every 
u G i?^(]R) let define the measure 

M"(A)=/ ulix)dx 

for every Borel set A C x T. The function ^'o can be regarded as a trans- 
portation map which transports the measure /x" into the measure n"" . The 
distance we shall introduce in Chapter 3 will be thus a sort of Wasserstein 
distance between measure [47]. 



Chapter 2 

Conservative multi-peakon 
solution 

"I was observing the motion of a boat which was rapidly drawn 
along a narrow channel by a pair of horses, when the boat suddenly 
stopped - not so the mass of water in the channel which it had put 
in motion; it accumulated round the prow of the vessel in a state 
of violent agitation then suddenly leaving it behind, rolled forward 
with great velocity, assuming the form of large solitary elevation, a 
rounded, smooth and well defined heap of water, which continued 
its course along the channel apparently without change of form or 
diminution of speed. I followed it on horseback and overtook it still 
rolling on at a rate of some eight or nine miles an hour, preserving 
its original figure some thirty feet long and a foot to a foot and a 
half in height, its height gradually diminished and after a chase of 
one or two miles I lost it in the windings of the channel. Such in 
the month of August 1834 was my first chance interview with that 
singular and beautiful phenomenon which I have called the Wave of 
Translation..." 

John Scott Russell, 1844 

It was Scott Russell [44] who introduces the concept of solitary waves to indicate 
no more than wave which propagate without change of form and have some 
localized shape (see also, [26, 27]). His experiment is described in fig. 2.1. For 
more than sixty years it was only a pure scientific curiosity, until Korteweg 
and de Vries [38] derived the equation for the propagation of waves in one 
direction on the surface of a shallow canal. The profile of the travelling wave 
solution with permanent shape found by Miura [40] is precisely the shape of 
the wave which Scott Russel observed in his experiments. The term soliton is 
substantially different from solitary wave, it was introduced in 1965 by Zabusky 
and Kruskal [50] to indicate waves that whenever collide each other they do not 
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/// 












Figure 2.1: Diagram of Scott Russell's experiment to generate a solitary wave 



break up and disperse, but remains almost identical to a solitary wave solution. 
In this chapter we investigate the shape of solitary waves for the Camassa-Holm 
equation (1.1). 



2.1 Multi-peakon solution in the real line 

In this section we shall construct a solution of the Camassa-Holm equation 
starting from an initial condition Uq of the form 

N 

The motivation of this choice is given by the shape of traveling wave solution 
(see [13] and [20, Example 5.2]). Looking for solution of the equation (1.1) in 
the traveling wave form u{t,x) = U{x — ct), with a function U that vanishes at 
infinity, the limit of k ^ leads to the function U = ce~l^~'^*l. This is not a 
solitary wave in the sense introduced by Scott Russell because of the presence of 
the cuspid at the position x = ct. However, the evolution of an initial data like 
uq remains of the same shape [13, 24]. It is a superposition of peaked solitary 
waves, which evolves as 

N 

«(i,a:) = ^ft(i)e-l--«^(*)l. 
i=i 

Hence we term peakon a peaked solitary wave to emphasize the soliton properties 
of such a function. 

As long as the classical solution of the problem 

( ^ 

(2.1) 

Pi = Pi^Pj sign(g'i - gijOe"!**"*^! 
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exists, the solution of this system gives the coefficients p{t) = (pi, . . . ,pn) and 
q(t) = (gi, . . . , qn) for the solution u{t, x) to the Camassa-Holm equation. Let 
us observe that the previous system can be viewed as an Hamiltonian system 
with Hamiltonian function i7(q,p) = i j Pi-Pj'^ ''^ ' • 

In [33] the authors prove the existence of a global multipeakon solution 
when strengths pi are positive for all i = 1 ... A'' and the convergence of the 
sequence of multipeakon solution. If uq is an initial data such that the dis- 
tribution Wo — woaa; is a positive Radon measure, there exists a sequence of 
multipeakons that converges in L°°{M., Hl^^{M.)). In this case the crucial fact is 
that no interaction between the peakons occurs, and then the gradient remains 
bounded. However, a general initial data contains both positive and negative 
peakons, as in the so called peakon-antipeakon interaction: one positive peakon 
with strength p, centered in —q, moves forward and one negative anti-peakon 
in q, with strength p moves backward (fig. 2.2). The evolution of the system 
produces the overlapping of the two peakons at finite time t = t, so that q ^ 0. 
The conservation of the energy E = H{q{t), p{t)) yields 



E 



lim p2(i_e-2|9l). 



(2.2) 



and then the quantity p blows up in finite time. At the point (r, 0) occurs 
thus a singularity for the solution u. To extend the solution also after the 

interaction time with a solution which conserves the energy E we can think that 
at the interaction point an antipeakon/peakon couple emerge, the first, negative, 
moving backward and the second, positive, moving forward with coefficients 

(— q, — p) and {q,p). According to the conservation of the energy, the choice of q 



and p must satisfy (2.2) as t 
the singularity at (r, 0) 



T+ . It yields a change of variables which resolves 



c 



(jj = arctan(p) 



with this choice, the Hamiltonian system leads to the ODE 



It 



(r)=( I 

2 



with 



1 _ e-Ccot=(a,) _ ^cot2(a;)e-^'^°t'('^) tan3(a;) 
sin2(a;)e-'^'^°*'(") 



and / is a Lipschitz vector field in a neighborhood of the point (-j, The 
solution {({t),uj{t)) of this problem provides then the unique couple {q{t),p{t)) 
which coincides with the classical solution of the Hamiltonian system for t < t 
and extends it for t > t. 

This example suggests the way to construct the multipeakon solution when- 
ever an interaction between peakons occurs (see also [10] for an "energetic" 
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t < T 
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Figure 2.2: Peakon-antipeakon interaction 



motivation). Suppose that two or more peakons with strengths p\,...,pk an- 
nihilate at the position q at time r and produce a blow up of the gradient ^i^- 
The conservation of the energy yields that there exists and is positive the limit 

/•«+(*) 

Br = lim / {t, x) dx 

where and ^+ are the smallest and the largest characteristic curve passing 
through the point (r, q). Assume that after the interaction two peakons appear 
with strengths pi, p2 and placed at the position qi, q2. Let consider the change 
of variables 



z = p2 + Pi w ^ 2 arctan(p2 ^ Pi) V 
then the system (2.1) turns out to be 

W = - \sm{w) cosh ( 2tan^^(^/2) ^ 

z'^ cos'^{w/2) - sin^(w;/2)]e~'-"''(™/2) 
isin(u;)sinh(2t5^ 



q2 + qi Q = {P2 - Pif {q2 ~ qi) 

) 



2^sinh L^^J(^^,) 



o>k+i 



{w/2)) + ^COSh (^2tan2^(«;/2) 



?7= z[l + e tan4»/2)] + 2 cosh 
- tan3(i 



/■ — p tan^(™/2) 

^ tan(tu/2) 



+2C 



2tan2(^i;/2 



i>/c+l 
5. 



sinh 



tan"'(u)/2) 
C 

— cosh I 



tan2(i«/2) 



tan(ti!/2) 
2 tan2^(u)/2) 



j>fe+l 



2.2 Approximation of the initial data 
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Pi = P,e [2 cosh ( ztanW) ) + tan(w;/2) sinh [ ^.J^^/^^ ) 



^2tan2(u)/2) 

+ E PiPjsign(gi-gj)e-l«'-«il 
j>fe+i 

-9i+'7/2 



+ 



j>fc+i 



cosh ( 2tan^'^(^/2) ) + tan(«;/2) sinh (atiJl^^ 



which is a system of ODE with locally Lipschitz continuous right hand side that 
can be extended smoothly also at the value w = w. The initial data become 

k 

z{t) = lim y]pi(i) w{t) = tt 77(t) = 2q ({t) = e^- 
1=1 

Piir) = lim Pi (i) qiir) = lim i = k + l,...,N 

Thus there exists a unique solution of such a system which provides a mul- 
tipeakon solution defined on some interval [r, r'[, up to the next interaction 
time. As wc will show in Corollary 3.1, since Camassa-Holm equation is time 
reversible, once we prove the uniqueness of the solution of a Cauchy problem, 
we have that maximal the number of peakons interaction is actually k = 2, one 
with positive strength the other with negative one. 

2.2 Approximation of the initial data 

In this section wc shall construct an approximation with initial data with a 
multipeakon function. Our aim is to approximate it with a sequence Ue which 
has an exponential decay at infinity uniformly w.r.t e. 



Lemma 2.1. Let f S Xa 

function g of the form 



Then for every e > there exists a multipeakon 

N 



g{x) = 5^fte-l--«*l 



=1 



such that 



\\f -g\\m{V) < e 

[g\x)+gl{x)]e''\-\dx<Co 



(2.3) 
(2.4) 



for some constant Co > which does not depend on e. 
Proof. Let p{x) G Cg" be a cut-off function such that 

• pix) > 

• p{x) = 1 for every |a;| < 1, p{x) = for every |a;| > 2 



J^p{x)dx = 1 
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and Pe{x) = -p{-) be a mollificrs sequence. Observe that for every e > 0, 
fe{x) = pe * f{x) is a smooth function which approximates the function / in 
ff^— norm 

\\f - fe\\m{M) < Cs (2.5) 
moreover it belongs to X^, indeed 



[[f^{x)+fel{x)]e»^^^dx < 
Jr 



s«l^l dx 



{f\y) + fl{y))Pe{x-y)dy 

< I [f\y) + fliv)] I Pe{x - 2/)e«l-l dxdy 

Jm. Jk. 

< I [f\y) + /J(2y)]Coe"l^l dy = CoC^'f < oo 
Jr 

and Co is a constant which does not depend on e. From the previous inequality 
we can assert that for every R > one has ||/£||i/i(R\[-H,iij) < CoCae~"^ 
uniformly in e > 0. We can choose thus big enough in order to have 

ll./ellHi(K\[--R.,«.]) < '^Z^- (2-6) 

In the space H^{[—Rs, R^]) we can approximate with a multipeakon func- 
tion. By using the identity 



-Jl-^]e-^^^=5o inV 



d_ 

2 V' dx 

the function can be rewritten in convolution form 

In the interval [—Re, Re] the previous integral can now be approximated with a 
Riemann sum 



N \ q^ = ^Re 



i=-N 



Pi= [ 

J Qi-l 



dy. 



Choosing N sufficiently large we obtain \\fe — fl||ffi([-_Re,i?e]) < ^- Together with 
(2.5) and (2.6) this last estimate yields the result. □ 

2.3 Periodic multi-peakon 

By a periodic peakon we mean a function of the form 

u{x)=vx{x-q), x(^) = E^"''""'- (2.7) 



2.3 Periodic multi-peakon 
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Observe that the periodic function x satisfies 

x{-x) = x{x) = x{x + l) xeR, 
, , e^ + e^-^ (2.8) 

e — 1 

For future use, we observe that for every periodic function u, the convolution 
P, defined in (1.2), takes the form 

P{x) = x{x - V) {u\v) + ^) dy (2.9) 

Wc begin this section by observing that also any periodic initial data can be 
approximated by periodic multi-peakons. 

Lemma 2.2. Let f e -ffper- Then for any s > there exists periodic multi- 
peakon g, of the form 

N N 

9ix) = E e-l^-"*-"! =J2pi X{x - qi) 

such that 

Proof. By taking a suitable mollification, we can approximate / with a periodic 
function / e C°°, so that 

<e/2. (2.10) 

Next, we observe that 

where 5o denotes the Dirac distribution concentrating a unit mass at the origin. 
We can thus write / as a convolution: 







2 

The above integral can now be approximated with a Ricmann sum 



i=l J{i-l)/N 



dy. 



Choosing sufBciently large we obtain ||/ — ffHiji^^ < ff/2. Together with (2.10) 
this yields the result. □ 
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Next, wc show how to construct a unique conservative solution, for multi- 
peakon initial data. As long as the locations qi of the peakons remain distinct, 
this can be obtained by solving the Hamiltonian system of O.D.E's (1.7). 

However, at a time r where two or more peakons interact, the correspond- 
ing strengths pi become unbounded. A suitable transformation of variables is 
needed, in order to resolve the singularity and uniquely extend the solution 
beyond the interaction time. 

Lemma 2.3. Letu be any periodic, multi-peakon iniiial data. Then the Cauchy 
problem (1.1), (l-H) has a global, conservative multi-peakon solution defined for 
all t & R. The set X of times where two or more peakons interact is at most 
countable. Moreover, for all t ^I, the energy conservation (1-13) holds. 

Proof. The solution can be uniquely constructed by solving the hamiltonian 

system (1.7), up to the first time r where two or more peakons interact. We 
now show that there exists a unique way to prolong the solution for i > r, in 
terms of two outgoing peakons. To fix the ideas, call 

q= lim qi{t) i=l,...,k, 

the place where the interaction occurs, and let pi{t), . . . ,Pk{t) be the strengths 
of the interacting peakons. Later in Section 3.5 we will show that only the case 
k = 2 can actually occur, but at this stage we need to consider the more general 
case. We observe that the strengths Pfe+i, . . . ,Pn of the peakons not involved 
in the interaction remain continuous at time r. Moreover, by (1.7) there exists 
the limit 



We can thus write 



P = lim_^pi(t) . 



N N 



u{t,x) = lim_^pi(i)e-l^-9'WI = pe-l^-«l+ ^ ^^(t) e-l^-«*(^)l . 

i=l i=k+l 

For t > T, we shall prolong the solution with two peakons emerging from the 
point q. The strength of these two peakons will be uniquely determined by the 
requirement of energy conservation (1.13). 

Call £,~{t), £,~^{t) respectively the position of the smallest and largest char- 
acteristic curves passing through the point {T,q), namely 

r(0 =min{e(t); ar)=q, i{s) = u{s,^{s)) for all s € [r - /i,T + /i]} , 
^+{t) = max|^(t); ^(r) = q, i{s) = u{s,^{s)) for all s e [r - /i, r + . 

(2.11) 

Moreover, define 

e(T,g)= lim / ul{t,x)dx. 
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The existence of this hmit follows from the balance law (1.5). This describes how 
much energy is concentrated at the interaction point. For t > t the solution 
will contain the peakons Pk+i, ■ ■ ■ ,Pn, located at gfc+i, . . . , Qn, together with 

the two outgoing peakons pi,p2, located at qi < q2- The behavior of pi,qi for 
i e {fc + 1, . . . , N} is still described by a system of O.D.E's as in (1.7). 

However, to describe the evolution of Pi,P2,qi, 52 one has to use a different 
set of variables, resolving the singularity occurring at (r, g). As f — > r+ we 
expect (1.8), (1.10) to hold. To devise a suitable set of rescaled variables, we 
observe that, by (1.3), 



along any characteristic curve t i-^ ^{t) emerging from the point q. Since u, P 
remain uniformly bounded, one has 



Ux{t,x)^- t>T, X G [qi{t) , q2{t)] . 

t — T 

The total amount of energy concentrated in the interval between the two peakons 
is given by 

r^'^ [u\t, X) + ul{t, x)] dx ^ ('^M_^\ ' ■{q2-qi) 
Jqi{t) V Q2-qi J 

(ps-pOll-e-l'^-^il)]^ 

q2-q\ 

« (P2 -Pl) (92 - gi) « (i{r,q) ■ 



The previous heuristic analysis suggests that, in order to resolve the singu- 
larities, we should work with the variables 

z=Pi+P2, = 2arctan(p2 -Pi) , 'n = Q2 + qi, C = {P2 - Pif {q2 - qi) , 

together with p^+i , . . . , pjv , gfe+i , • • • , gjv- 

To simplify the following calculations we here assume Q < qi < q2 < qk+i < 
... < qN < which is not restrictive. 

Let X defined in (2.7) and x(x) = ~^ — , x € (0, 1). From the original 
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system of equations (1.7) it follows 



i = cosh ( - cot^ ti;/2 j 2:cos^ pjx{qj — r]/2) 



N 



sinh I ^ cot^ ui/2 j tanii;/2 ^ pjx{qj — ri/2) 



j=k+i 

w = {z^ cos^ w/2 - sin^ w/2) x (Ccot^ w/2) 

+2 cosh [ - cot^ w/2 j z ^ pjX {qj - ?7/2) 

+2sinh ( - cot^«;/2 j sinw pjxiqj — rj/^) 



j=k+i 



r] = z 



[x(0) + X (C cot^ w;/2)] + cosh ( | cot^ «;/2 ) ^ " r?/2) 



3=k+l 



C = [x(0) - X (C cot^ w/2)] C tan^ «,/2 + x w/2^ z'^Q cot to/2 
— sinh f - cot^ w/2 j tan^ t(;/2 ^ pjx(gj— ry/2) 

cosh I - cot^'u;/2 I 2cotw/2 pjxiqj — t]/'^) 



+2Ccot'u;/2 



■ sinh cot^ w/2^ Yl PjX ilj - ri/2) 



Pi = Pi 



c 



cosh - cot w/2 I zx (qi — ri/2) + sinh I - cot w/2 I tan w/2x (<7i — ri/2) 



c 



AT 



+Pi XI Pi sign(gj - qj)x {\qi-qj\) 

j=k+i 

c 



= cosh I - cot w/2 I zx {qi — r]/2) + sinh I - cot w/2 I tan w/2x (ft — ??/2) 



c 



AT 

+ Yl Pjx{\qi-qj\) 
i=fe+i 

with initial data 

z{t)=p, w{t) =Tr, r]{T) = 2q, Cir) = e^r,q) , 

Piir) = lim pi{t), qiiT) = \im qi{t) i = k + l,...,N . 

For the above system of O.D.E's, a direct inspection reveals that the right 
hand side can be extended by c;ontinuity also at the value w = tt, because all 
singularities are removable. This continuous extension is actually smooth, in a 
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neighborhood of the initial data. Therefore, our Cauchy problem has a unique 
local solution. This provides a multi-peakon solution defined on some interval 
of the form [r, t'[ , up to the next interaction time. 

The case where two or more groiips of pcakons interact exactly at the same 
time r, but at different locations within the interval [0,1], can be treated in 
exactly the same way. Since the total number of peakons (on a unit interval 
in the .t- variable) does not increase, it is clear that the number of interaction 
times is at most countable. The solution can thus be extended to all times t > 0, 
conserving its total energy. □ 
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Chapter 3 

Distance defined by optimal 
transportation problem 



3.1 A distance functional in the spatially peri- 
odic case 

In this section wc shall construct a functional J(u, v) which controls the distance 
between two solutions of the equation (1.1). All functions and measures on M 
are assumed to be periodic with period 1. Let T be the unit circle, so that 
T — [0, 27r] with the cndpoints and 2tt identified. The distance \d — 
between two points 9,9' G T is defined as the smaller between the lengths of 
the two arcs connecting 6 with 9' (one clockwise, the other counterclockwise). 
We now consider the product space 



where a Ab = min{a, b}. Let Ai{X) be the space of all Radon measures on X 
which arc 1-pcriodic w.r.t. the .r-variable. To each l-periodic function u € -ffper 
we now associate the positive measure cr" G ^A{X) defined as 



for every Borel set ^ C x T . Notice that the total mass of a^"''*-* over one 
period is 



X = M X K X T 



with distance 




(3.1) 




(3.2) 
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On this family of positive, 1-pcriodic Radon mcasnrcs, wc now introduce 
a kind of Kantorovich distance, related to an optimal transportation problem. 
Given the two measures tr" and <t", their distance J{u,u) is defined as follows. 

Call the family of all strictly increasing absolutely continuous maps ip : 
R I— > M which have an absolutely continuous inverse and satisfy 

ipix + n) = n + tp{x) for every n e Z . (3.3) 

For a given € JF, we define the 1-periodic, measurable functions 0i , 02 : M i-^ 
[0, 1] by setting 



(3.4) 



</>i(x) = sup 1^ G [0, 1] ; e-(l + ul{x)) < (l + uli^ix))) ^'{x) 
ct>2{ip{x)) = sup |e e [0, 1] ; 1 + ul{x) >e-(l + ul{ip{x))) V''(a;)| . 

Observe that the above definitions imply max |(/)i(a;), (^2(2^)} = 1 together with 
Mx)(l + ul{x))=^2{ip{x)) (l+ul{i,{x)))i,'{x) (3.5) 
for a.e. a; e K. We now define 

■(l)i{x) (1 + u1{x)) dx 

+ f (1 + ul{x)) - (1 + ul{ij{x))) ij'{x) dx . 
Jo 

(3.6) 

Of course, the integral is always computed over one period. Observe that the 
map X 'tpix) can be regarded as a transportation plan, in order to transport 
the measure cr" onto the measure cr". Since these two positive measures need 
not have the same total mass, we allow the presence of some excess mass, not 
transferred from one place to the other. The penalty for this excess mass is 
given by the second integral in (3.6). The factor < 1 in the first integral 
indicates the percentage of the mass which is actually transported. Integrating 
(3.5) over one period, we find 

/ Mx){^ + ul{x))dx= [ My){^ + uliy)) dy . 

Jo Jo 
We can thus transport the measure (pi cr" onto ^2 by a map 

^ : (x, u{x) arctanMa;(a;)) i— > (j/, u{y), a.rcta.nux{y)) , 

where y — ip{x). The associated cost is given by the first integral in (3.6). 
Notice that in this case the measure 02 is obtained as the push-forward 
of the measure 0i cr". We recall that the push-forward of a measure cr by a 
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mapping ^ is defined as {^^(7){A) = ^(^)) for every measurable set A. 
Here = {z ; -^{z) e A}. 

Our distance functional J is now obtained by optimizing over all transporta- 
tion plans, namely 

J{u,u) = inf J^(u,{t). (3.7) 

To check that (3.7) actually defines a distance, let u,v,w & -ff^(M) be given. 

1. Choosing ^(x) = x, so that (t>i{x) = 4>2{x) = 1, we immediately see that 
J{u,u) = 0. Moreover, if J{u,u) = 0, then by the definition of d'^ we have 
u = u. 

2. Given ij) & define tp = tp~^, so that = 02, ^2 = <^i- This yields 

J'^{u,u) = J'l'{u,u). 

Hence J{u,u) = J{u,u). 

3. Finally, to prove the triangle inequality, let -(/^^ , -0" : R i— !■ R be two increas- 
ing diffcomorphisms satisfying (3.3), and let (j>\, (j)^, (j>\ : M [0,1] be the 
corresponding functions, defined as in (3.4). We now consider the composition 
tp = ipi otp^ and define the functions 0i, (/>2 according to (3.4). Observing that 

Mi'i^)) = <^2(v"(/(a;))) > <^^(V'(x)) •4(v«(/(x))) , 

and recalling that the distance d'^ at (3.1) is always < 1, we conclude 

(v,w) . 

This implies the triangle inequality J{u, v) + J{v, w) > J{u, w). □ 
In the remainder of this section we study the relations between our distance 
functional J and the distances determined by various norms. 

Lemma 3.1. For any u,v G -ffper one has 

with a constant C uniformly valid on hounded subsets of H^^^. 
Proof. We shall use the elementary bound 

I arctana — arctan6| • < 47r(|a| + |6|) |a — fe| , 
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Figure 3.1: Infinitesimal area between u and v 



valid for all a, 6 e K. In connection with the identity mapping ^|J{x) = x we now 
compute 



J'^{u,v) < / \\u{x)~v{x) 
Jo ^ 

Jo 



+ 2 1 arctan-Ua; — arctan 
dx 



dx 



< \\u- v\\loo \\1 + uJ\li + {Sn + 1) / \ux +Vx\\ux - Vx\dx 

,, . Jo 

< (Stt + 3) (1 + \\u\\h^ + \\v\\h^) ■ \\u - v\\h^ , 

proving the second inequality in (3.8). 

To achieve the first inequality, choose any ip G J^. For x G [0, 1], call 7^ the 
segment joining the point = [x,u{x)) with = (^ip{x),v{tlj{x))) . Clearly, 
the union of all these segments covers the region between the graphs of u and 
V. Moving the base point from x to a; + dx, the corresponding segments sweep 
an infinitesimal area dA estimated by (fig. 3.1) 



\dA\ < |p--g-|.(|dP"| + |ciQ-|) 

< {\x - V(x)|' + \u{x) - v{^P{x))\'y^'\{l + ulfl^ + (1 + vlfl^\x) 



dx . 



Integrating over one period we obtain 

\u{x) - v{x)\ dx < j - ■0(t)(a;)| + |u(x) - u(-0(t)(a;))|^- 

■ [(1 + ulix)) ^'^ + (1 + vl{^{x))) V'w {x) 

< {2+\\u\\H^ + \\v\\H^)-[j^iu,v) + J^-\v,u) 

<C-J{u,v), 
completing the proof of (3.8). 



dx 



□ 
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Lemma 3.2. Let (u„)„>i he a Cauchy sequence for the distance J, uniformly 

bounded in the H^^^ norm. Then 

(i) There exists a limit function u € -ffpor such that Un u in L°° and the 
sequence of derivatives Un^x converges to Ux in ipg^; for 1 < p <2. 

(ii) Let Un be the absolutely continuous measure having density ^ with re- 
spect to Lebesgue measure. Then one has the weak convergence jin !->■, 
for some measure fj, whose absolutely continuous part has density . 

Proof. 1. By Lemma 3.1 we already know the convergence u„ — > u, for some 
limit functfon u G ^^per ■ By a Sobolev embedding theorem, all functfons m„,u 
are uniformly Holder continuous. This implies ||uri — — > 0. To establish 
the convergence of derivatives, we first show that the sequence of functions 

Vn = exp{2iarctanu„_^} 

is compact in L^^^ . 

Indeed, fix e > 0. Then there exists N such that J{um, u„) < e for m,n > N . 
We can now approximate M^r in H^^^ with a piecewise affine function u^v such 
that J{un,Un) < s. By assumption, choosing suitable transport maps Vn we 
obtain 



Jo 



exp {2iarctan'u„^a,(a;)} — exp {2iarctanujv,x(V'n(3;))} dx < 2J{un,iiN) 

< Ae 



for all n > A^. We now observe that all functions 

X I— > exp {2iarctanujv,x(V'n(2;))} 

are uniformly bounded, piecewise constant with the same number of jumps: 

namely, the number of subintervals on which un is affine. The set of all such 
functions is compact in L^^^. This argument shows that the sequence v„ = 
exp{2i arctan w„^a;} eventually remains in an £-neighborhood of a compact subset 

of ipci-. Since e > can be taken arbitrarily small, by possibly choosing a 
subsequence we obtain the strong convergence Vn v for some v G L^^^ . 

2. From the uniform bounds and the L^ convergence of the hmctions v„, 
we now derive the L^ convergence of the derivatives. For a given e > 0, define 

M = sup||«„||h-i , An = \x e[0,l]; \un,x{x)\ > M/e\ . 

^ per (_ ' ' ) 

The above definitions imply 

meas(^„) < £^ 
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We now have 

\\Um,x - Un,x\\LP < 




[0,1]\(A„UA„) 



\'Pdx 



l/p 



i/p 



h < 







= h+h. 






(2-p)/2p 


< 


/ 1-dx 










< e(2-p)/p . 2M. 





/ {\Um,x\ + Wn,x\T dx 

.Ja^ua„ 



1/2 



(3.8) 



Next, choosing a constant such that 



_2i arctan a Jli arctan 6 1 



>Ce|a-6| 



we obtain 



arctan g2i arctan ii^^j; I P 



whenever |a|, |6| < Mje , 
i/p 



(3.9) 



Taking e > small, wc can make the right hand side of (3.8) as small as we 
like. On the other hand, choosing a subsequence such that = e^*'^'''^**^'^"''.* 
converges in Lpgj., the right hand side of (3.9) approaches zero. Hence, for this 
subsequence, 

limsup \Um,x - Wn.xlUg,, = . 



m.n—^oo 



Since Un ^ u uniformly, in this case we must have 

\\Un,x - UxWl^,,^ . 



(3.10) 



We now observe that from any subsequence we can extract a further subsequence 
for which (3.10) holds. Therefore, the whole sequence (un,x)n>i converges to 
Ux in L^er • 

3. To establish (ii), we consider the sequence of measures having density 1+m^ ^ 
w.r.t. Lebcsguc measure. This sequence converges weakly, because our distance 
functional is stronger than the Kantorovich-Waserstein metric which induces 
the topology of weak convergence on spaces of measures. Therefore, /i„ — ^ /x for 
some positive measure /i. 

Since the sequence 1+Un,x converges to 1+Ux in Lp^r, by possibly choosing a 
subsequence we achieve the pointwise convergence Un^x{x) Ux{x), for a.e. x G 
[0, 1]. For any e > 0, by Egorov's theorem we have the uniform convergence 
Un,x{x) Ux{x) for all x G [0, 1] \ T^, for some set with meas(y£) < e. Since 
£ > can be taken arbitrarily small, this shows that the absolutely continuous 
part of the measure n has density + w.r.t. Lebesgue measure. □ 
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3.1.1 Continuity in time of the distance functional 

Here and in the next section we examine how the distance functional J(-,-) 
evolves in time, in connection with multi-peakon solutions of the Camassa- 
Holm equation (1.1). We first provide estimates valid on a time interval where 
no peakon interactions occur. Then we show that the distance functional is 
continuous across times of interaction. Since the number of peakons is locally 
finite, this will suffice to derive the basic estimates (1.15)-(1.16), in the case of 
multi-peakon solutions. 

Lemma 3.3. Let t uit) G -ffpor rnulti-peakon solution of (1-1). Assume 
that no peakon interactions occur within the interval [0,t]. Then 

J{u{s),u{s')) <C -{s- s'\, for all s,s' &[Q,t], (3.11) 

for some constant C, uniformly valid as u ranges on hounded subsets of H^^t ■ 



Proof. Assume Q < s < s' < t. By the assumptions, the solution u = u(t,x) 
remains uniformly Lipschitz continuous on the time interval [0, r]. Therefore, 
for each s e [0, r] and a; e M, the Cauchy problem 

determines a unique characteristic curve t ^(f; s, .t) passing through the point 
{s,x). Given s' G [0,t], we can thus define a transportation plan by setting 

'il,{x)=^{s';s,x). (3.13) 

Of course, moving mass along the characteristics is the most natural thing to 
do. We then choose , (/)2 to be as large as possible, according to (3.4). Namely: 

(x) = sup 1^ e [0, 1] ; 6-{l + ul{s, x)) < i>'{x) -(l + ul {s' , | , 

Mx) = sup 1^ e [0, 1] ; ^ • (l + ul{s', V(a;))) i>'{x) < 1 + ul{s, a;)| . 
The cost of this plan is bounded by 

J''' {u{s) , u{s')) < [ \\x-^{s';s,x)\ + \u{s,x)-u{s',^{s';s,x))\ 

+ |2arctan'Ux(s,a:;) — 2arctan'U2;(s',^(s';s,a;))|^ Ml + u'^{s,x)) dx 

+ [ (1 -0i(a;)) {l + ul{s,x))dx 

+ j (1-02(V(^))) (l + «^(s',V(^)))^'(x)rfa;. 

(3.14) 
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To estimate the right hand side of (3.14), we first observe that, for all u G -ffpen 



li^cc ^ J \u{x)\dx + J \ux{x)\dx 

<ML^ + \\u4i^^<2\\u\\H^^^ = 2iE-y/'. 



Using (3.15) in (3.12) we obtain 

\^{s)-^{s')\<2{Ey/' 



s-s- . 



Next, from the definition of the source term P at (1.2) it follows 



-\x\ 



< \\u\\ 



Similarly, 
Using (3.18) we obtain 

*(s', as')) - u{s, as)) < [ 

J s 



LH[0,1]) 

1 = 



per 



(3.15) 

(3.16) 

(3.17) 
(3.18) 



d 



<t, m) 



Px{t, m) 



dt 

dt < E""- Is' - s\ . 



dt 



(3.19) 



Concerning the term involving arctangents, recalling (1.3) we obtain 
d 



dt 



2 arctanwa;(f, ^(t, a;)) 



2 




~ 1 + ul 


2 



The bounds (3.17) and (3.19) thus yield 

2arctanUa;(s',^(s')) — 2arctanua;(s, ^(s)) 

< (2\\u\\l^ + l + 2\\PU^)-\s'-s\ (3-20) 

< (I0i;" + 1) • \s' -s\. 

This already provides a bound on the first integral on the right hand side of 
(4.4). 

Next, call /i, /2 the last two integrals on the right hand side of (4.4). Notice 
that 



Ii+l2 = {l+ul{s,v))-^y{s';s,y){l+ul{s',as';s,y))) 



dy 



Indeed, Ji + I2 measures the difference between the measure (l + w^(s',2/)) dy 
and the push- forward of the measure (l + u'^{s,x))dx through the mapping 
X ^(s';s,a:). 
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Since the push-forward of the measure u| dy satisfies the linear conservation 



law 

wt + {uw)a, = , 
comparing (3.21) with (1.5) we deduce 



(3.21) 



Jo 



ul{s,y)-^y{s,y)ul{s',^{s';s,y)) 

J s JO 

< 2 I (\\u\\l^ + \\P\\l^) Wu^hi dt < 2 UE^^ + E^) E^ ■\s'-s 



dy < I [ 2\{u^ - P)u^\dxdt 

s Jq 



because of (3.15), (3.17) and (1.13). Finally, we need to estimate the remaining 
terms, describing by how much the Lebesgue measure fails to be conserved by 
the transformation x ^{s' ; s, x). Observing that 



d_ 

dt 



^y{t,y) =u^{t, ^it,y))^y{t,y), 



^j/(0,y) = l, 



we find 



\'^-^v{s';s,y)\dy< J Q^[^y{t,s,y)] 
< j / iv{t;s,y)\u^{t,£,{t\s,y))\dydt. 

Js Jo 



dydt 



(3.22) 



To estimate the right hand side of (3.22), we use the decomposition [0,1] = 
YUY'UY", where 

Y={y; s, y) e [(1/2) , 2] for aU t e [s, s'] } , 

Y' = {y-, iy{t; s, y) < 1/2 for some t G [s, s'] | , 

y" = 1 2/ ; ^y(t; s,y) >2 for some t G [s, s'] | . 
Integrating over Y one finds 

/ / ^yit;s,y)\u^{t,^{t;s,y))\dydt < 2 [ \\u^{t)\\ dt < 2^" • |s' 

Js Jy Js 

Next, if y G y we define 

T{y) =M {t> s; ^y{t;s,y) < 1/2} . 
Observe that y gY' implies 



Si . 



I Ux{t, ^{t;s,y)) 

J s 



dt > ln2. 
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Therefore 



dt 



dy 



< 



Ux{t, i{t;s,y)) 
^ [ [ \u^(t,x)\dxdt < 4:E^-\s'-s\ 



The estimate for the integral over Y" is entirely analogous, Indeed, the push- 
forward of the Lcbesgue measure along characteristic curves from t — stot = s' 
satisfies exactly the same type of estimates as the pull-back of the Lebesgue 
measure from t = s' to t = s. All together, these three estimates imply 



X 



YUY'UY" 



\^-^v{s';s,y)\dy< lOE^- \s' - s\ 



(3.23) 



Putting together the estimates (3.16), (3.19), (3.20), (3.13) and (3.23), the 

distance in (4.4) can be estimated by 

J"^ (w(s) , u{s')) < [2(1 + £")+£" + (10 i;" + 1) + 10 (E'^y + 10 E^] • |s' - s| . 
This establishes (3.11). □ 

According to Lemma 3.3, as long as no peakon interactions occur, the map 
1 1—!- u(t) remains uniformly Lipschitz continuous w.r.t. our distance functional, 
with a Lipschitz constant that depends only on the total energy £■". Since in- 
teractions can occur only at isolated times, to obtain a global Lipschitz estimate 
it suffices to show that trajectories are continuous (w.r.t. the distance J) also 
at interaction times. 

Lemma 3.4. Assume that the multi-peakon solution u{-) contains two or more 
peakons which interact at a time r. Then 

lim J(u(t - h), u(t + h)) = . 

Proof. To fix the ideas, call x = q the place where the interaction occurs, and 
let pi, ... ,pk be the strengths of the peakons that interact at time r. We here 
assume that < g < 1. The case where two or more groups of peakons interact 
exactly at the same time r, within the interval [0, 1], can be treated similarly. 

For \t — T\ < h, call ^~{t), £,^{t) respectively the position of the smallest and 
largest characteristic curves passing through the point (r, ^), as in (2.11). We 
observe that u is Lipschitz continuous in a neighborhood of each point (t, x) , 
with X q. Hence, for x G [0, 1] \ {q} there exists a unique characteristic curve 
ti-^^{t; T, x) passing through x at time r. For a fixed h > 0, the transport map 
ip is defined as follows. Consider the intervals I-h = [C~('''~ ^)i h)] and 

Ih = [$~(t + h), ^^{t + h)] . On the complement [0, 1] \ I-h we define 



V'(€(t - h; T, x)) = ^(r + h; r, x) , 
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so that transport is performed along eharaeteristic curves. It now remains to 
extends ^ as a map from I-h onto Ih- Toward this goal, we recall that our 
construction of multi-peakon solutions in Section 2.1 was specifically designed 
in order to achieve the identity 

e(^^q) = ^lim^ y u^{T — h,x)dx = ^im^ J ul{T + h,x)dx. (3.24) 
For ft. > we introduce the quantities 

E{-h)= {l+ul{T-h, x))dx, E{h)= I {1 + uI{t + h, x)) dx , 

Ji-h Jih 

e{h) = 2mm{E{-h), E{h)} - max{E{-h), E{h)} . 

Notice that (3.24) implies e{h) > and 

E{-h) e(r,q) , E{h) e(^r,q) , e{h) e^r,g) , (3.25) 
as /i — > 0+. Consider the point x* = x*{h) inside the interval 
I.h=[r{r-h), eiT-h)], 

implicitly defined by 

/ [l + u1.{x)) dx = e{h) . 

h-{T-h) 

For X e [^~(t — h), a;*] we define i}{x) as the unique point such that 

/■l/'(x) rx 

I {1 + uI{t + h, x)) dx = (1 + uI{t - h, x)) dx . (3.26) 

h-(T+h) h-{r-h) 

We then extend ij} as an affine map from \x*, ^"""(r — /i)] onto [^'(a;*), ^"""(r+Zi)] , 
namely 

i){9-£,+{T-h) + {i-e)-x*^ =e-^+{T+h) + {i-e)-i){x*) e& [o,i]. 

Finally, we prolong ip to the whole real line according to (3.3). 

As usual, the 1-periodic functions (f>i , 4>2 are then chosen to be as large as 
possible, according to (3.4). As /i — > 0+, we claim that the following quantity 
approaches zero: 



J'l'{u{T-h), u{T + h)) 

[tp{x), u{t + h, tp{x)), 2 arctan{(a,(r + h, ip^x))) 
■Mx) + uI{t - h, x)) dx ' (3-27) 

+ / {\- ^x{x)){\^uI{t -h,x))dx 

+ y (1 - 4>i{i'{x))) (l + (r + ft, i,{x))') V'(x) dx . 
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It is clear that the restriction of all the above integrals to the complement 
[0, 1] \ I-h approaches zero as /i ^ 0. We now prove that their restriction to 
I-h also vanishes in the limit. As /i — > 0+, for x G I-h we have 

(j^x, u{t — h,x), 2 arctanMa:(T — h, x)) , 

(^p{x), u{t + h,ijj{x)), 2ai,Tctai,nUx{T + h,'ijj{x)))^ 0, 

because all points approach the same limit (g, u{T,q), tt). The first integral in 
(3.27) thus approaches zero as /i ^ 0+. 

Concerning the last two integrals, by (3.26) it follows 

Mx) = 02 (^-(2;)) = 1 for all X G [C"(t - h), x*] . 

Moreover, our choice of x* implies 

/ {l + uI{t - h, x)) dx + [1+uI{t + h, ■ilj{x)))ilj'{x)dx 

< 2m&x{E{-h), E{h)} - 2mm{E{-h), E{h)} . 

(3.28) 

By (3.25), as /i — > 0+ the right hand side of (3.28) approaches zero. Hence the 
same holds for the last two integrals in (3.27). This completes the proof of the 
lemma. □ 



3.1.2 Continuity w.r.t. the initial data 

We now consider two distinct solutions and study how the distance J{u{t) v{t)) 
evolves in time. To fix the ideas, let 1 1— > u{t) and 1 1— > v{t) be two multi-peakon 
solutions of (1.1), and assume that no interaction occurs within a given time 
interval [0,T]. In this case, the functions u,v remain Lipschitz continuous. We 
can thus define the characteristic curves t ^ ^{t,y) and t i— > C{t,y) as the 
solutions to the Cauchy problems 

i = u{t,o, m = y, 
C = v{t,o, m = y, 

respectively. Let now ■i/'(o) G be any transportation plan at time t = 0. For 
each t e [0, T] we can define a transportation plan ■)/'(*) G by setting 

^{tMt,y)) =C{t, V(o)(2/)) 

The corresponding functions (j)^^ , (j)2^ are then defined according to definitions 
in(3.4), namely 

</.f^(x)=sup|0G[O,l]; 0-(l+ul{t,x)) < (l + vl{t,i;^t-,{x)))i;[,){x)\ , 
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(a;)=sup [0,1]; l + ul{t,x) > ^ • (l + V(t)(a;))) V(t)(a;)| • 



If initially the point y is mapped into y = ip(^o){y), then at a later time t > the 
point ^{t, y) along the u-characteristic starting from y is sent to the point C,{t, y) 
along the ^-characteristic starting from y. We thus transport mass from the 

point (^^{t,y) , u(t,£_{t,y)) , 2 arctanuj; (t, ^(t, j/))^ to the corresponding point 

(C{t,y), v{t,C{t,y)), 2arctani;,(t,C(i,y))) , 

In the following, our main goal is to provide an upper bound on the time 
derivative of the function 



J'*^*\u{t), v{t)) = J d'^(^{x, u{t,x), 2arctanu^(t,a;)) , 

{il>(t){x), v{t, '0(t)(a:)), 2arctanu^(f, '0(t)(a;)))) • (t>i\x) {l + ul{t,x)) dx 
+ [ {l-4\x)){l + ul{t,x))dx 

+ I {l + vl{t,^|;^^^){x)))i^[^■^{x)dx. 



(3.29) 

Differentiating the right hand side of (3.29) one obtains several terms, due to 

• Changes in the distance d'^ between the points u, 2arctan'Ua;) and 
(C, V, 2arctanUj;). 

• Changes in the base measures (1 + u^) dx and (1 + u^) dx. 

Throughout the following, by 0{1) we denote a quantity which remains uni- 
formly bounded as u, v range in bounded subsets of -ffpgr • Using the elementary 
estimate 

I'" — 'f^l < (l + \u\ + \v\) min {|u — v\, l} , 
we begin by deriving the bound 

h = -\x-^^t){x)\-cl>^i\x){l+ulit,x)) dx 

J u{t,x)-v{t,ip(t){x)) ■ (l)f'{x){l + ul{t,x)) dx 

(l + ||«(i)|Loo + II^^Wll^oo) • J^W(«(t), V{t)) 

0(1)- J'^'w {u{t),v{t)). 



< 
< 



Here and in the sequel, the time derivative is computed along characteristics. 
Next, recalling the basic equation (1.1), we consider 



- '"I 



u{t,x) - v{t, ip(t){x)) 



< 



■4>t\x) {l + ul{t,x)) dx 
■ (l + w^(t, x)) dx . 
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In the spatially periodic case, by (1.2) and (2.8) we can write the source terms 



pu pv „„ 



1 r 

= 2j ^ x'(a; - y) ■ 



+ 



x'(V'(i)(a;) - y) ■ v^{t,y) dy 

I x'{'^(t){x)-i^{t){y)) 

Jx-l 



dy, 



where, according to (2.8), 



e- 1 



< a; < 1 . 



x'ix) = x'ix + i) xeR. 



In the next computation, we use the estimate 

{l + uliy))-(l + vl{^{y))y{y) 



dy = 0{l)-J'^{u,v). 



which holds because of the last two terms in the definition (3.6). Observing that 
x' is Lipschitz continuous on the open interval ]0, 1[ , we now compute (omitting 
explicit references to the time t) 

P:{x) - P:{i^{x)) \<\l' \x'{x - y) ■ u\y) - x'(V(a;) - y) ■ v\y)\ dy 



+0{l) ■ x-y- {i,{x) - V(2/)) 



+ \ I x!{x-y)Ul{y)-vl{'4){y))xl}'{y))dy 

* J x — l 

= 0{l)-{\^x-^l;{x)\+\\u^-v^L.) 
+0{l)■i\x-^l,{x)\+ Jy- 



V'(y)| 



tp'iy) dy 



+0{1)- (j^{u,v)+ / x\x-y)-[^'{y)-l]dy 
0(1) • |x- V(a;)| +0(1) • J'l'{u,v) 



+0{1) 



Hx)\ 



x-l 



X"{x - y) ■ [xl){y) - y] dy 



(3.30) 



= 0(1)- |a;-^(x)| +0(1)- J*(u,t;). 

Integrating over one period we conclude 

h = 0{l)-J^{u{t), v{t)). 
For future use, we observe that a computation entirely similar to (3.30) yields 

P"(a;)-P''(V'(a;)) = 0(1) • |a; - V(a;)| + 0(1) • J'^(u, t;) . (3.31) 



3.1 A distance functional in the spatially periodic case 
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Next, we look at the term 



dt 



2 arctan {t. x) — 2 arctan Vx {t, ijj(t) (x)) ■ (f)i^ {x) (l + u^. {t, x)) dx . 



Along a characteristic, according to (1.3) one has 

|2arctann.(t,^(t)) = ^ 
Call 0^ = 2 arctan lia; , = 2 arctan ij^; , so that 



1 2 6'" 

1+^ = ^°^ T' 



Ux 1 . , 

TT = Sin ( 

1 + ul 2 



1 + ul 



■ sm 



We now have 







(l + ulix)) 



uUx) 



1 + m2(x) 1 + 7j2(V;(.x)) 



dx 



JO, 



. ^e^{x) . ^e-(^p{x)) 

sin — sin — — - 



dx 



(3.32) 



< j {l + ul{x))- 9''{x)-e'"{^{x)) dx 
= 0{1) ■ J{u, v) . 
Next, using (3.31) we compute 

'-{x)-P"ix) v^{^p{x))-P^{ij{x)) 



I {^ + ul{x)) 
Jo 

<-f 

Jo 



u^{x) 



l + ul{x) ^ l + vl{i^{x)) 

v'^{il){x)) dx+ [ P"(a;) -P''(V'(x)) 
Jo 



dx 
dx 



+0{l) ■ f 
Jo 

0{l)-J^{u,v), 



l + ul{x) l + vl{i,{x)) 



(l + u^{x)) dx 



where the last term was estimated by observing that 

'< 12 arctan — 2 arctan 



1 + 1 + ^2 

Putting together all previous estimates wo conclude 

li+h+h=0{l)- J'^{u,v). 



To complete the analysis, we have to consider the terms due to the change 
in base measures. From (1.5) it follows that the production of new mass in the 
base measures is described by the balance laws 

r (1 + ui)t + [u{i + uD] ^ = [2u' + 1 - 2P^]ux = r , 

\ (1 + t;^), + [^,(1 + vl)\^ = [2^;2 + 1 - 2P^]vx = F • 
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This leads us to consider two further integrals h,!^ : 
h 



2 arctan?i2;(x)) , (^ip{x), !;(■!/; (x)), 2 a.rcta.nvx{'4'{^))) 
■\2u'^{x) + 1 - 2P"(x)| \u^{x)\ dx 



0{l)- J'l'iu^v). 



•(l + u^{x)) dx 



h = 



< 2 



/ [2u^{x) + 1 - 2P"(x)]m,(x) - [2i;2(^(x)) + 1 - 2P''(V(x))li;4^(x)) ^'{x) 

Jo L J 

jjl' l^\u\x)-v\'il;{x))\ + \p^ix)-Pm;{x))\^\u4x)\dx 

I 2t;2(^(x)) + 1 - 2P\^l^{x)) ■ Ux{x) - v,{^lj{x))^l;'{x) 
Jo 



dx 



Using (3.31) we easily obtain 



I', < 







u 



(x) ^ v"^ {i;{x)) + P'^ix)- P''{i^{x)) \{l + ul{x))dx 



0{1) ■ J'l'{u,v). 



On the other hand, recalling (3.32) and using the change of variable y = tp{x), 
x = tlJ~^{y), we find 



= 0{1) ■ f 
Jo 

= 0{1)- f 
Jo 



Ux{x) - Vx[ip{x))ip'{x) dx 



(l + u^{x)) dx 
l+ul{x) 



l + ul{x) l + vl{i,{x)) 

+0{l)- f v.4iP(x))^'{x)' 
Jo 

< 0{1) ■ J^{u, v) + J^'\{l + ulix)) - (1 + vli^ix))) i,'{x] 
= 0{1)-J{u,v). ° 



{l+vl{i,{x)))i,'{x) 



dx 



d_ 
It 



All together, the previous estimates show that 

J'^w(u(t), v{t)) < h+h + h+h + I'^ + I'^ = 0(1)- J'^«(«(i), v{t)) , (3.33) 



where 0(1) denotes a quantity which remains uniformly bounded as u,v range 
on bounded sets of -tfpgr • As an immediate consequence we obtain 



3.2 A priori bounds 
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Lemma 3.5. Let t u(t), t v(t) be two conservative, spatially periodic 
multipeakon solutions, as in Lemma 2. Then there exists a constant k, depending 
only on max{||u||Hi^^ , such that 

J{u{t) , v{t)) < e"'*-"! • J{u{s) , v{s)) s,tGR. (3.34) 

Proof. For t > s the estimate (3.34) follows from (3.33), taking the infimum 
among all transportation plans -0(5) at time s. The case i < s is obtained 
simply by observing that the Camassa-Holm equations are time-reversible. □ 

3.2 A priori bounds 

In [19] the author discusses the finite propagation speed property for the Camassa- 
Holm equation. Due to the nonlocal nature of the equation 1.1, it is not a priori 
clear that the evolution of an initial data with compact support will remains 
with compact support. On the contrary, [19] prove that the finite propagation 
speed property is valid only for the function and not for u. In this 

section we start from this result, and we want to estabilish what is the "right" 
decay at the infinity of solutions to 1.1. For this purpose, we introduce the 
following functional space. Let a e]0, 1[, then we set 

Xa={u€ H^{R) s.t. C"'" = / [u^{x) + ul{x)\ 6"'"=' dx < +oo}. (3.35) 

The present section is devoted to the study of some \iscful properties of the 
functions u £ X^. We start recalling an estimate for the norm of the 
fl'i(M) functions. We have 

< l|/|ll^i(M). (3.36) 

This estimate give us a bound on the norm of the conservative solution u 
of (3), in fact the conservation of the energy yields 

\\u{t)\\L^ < \\u{t)\\H^R) = ^ for every t > 0. (3.37) 

Let us consider now the behaviour of the functions u € Xa as |a:| goes to infinity. 
If we denote with C"'" the constant /jj(u^ + M^)e"'^' dx, the following holds 

sup M2(a;)e«l^l < 2C"'" . (3.38) 

Indeed, the function 

/(a;) =M(t,a;)e^l^l 

belongs to iJ^(M), moreover 

fx = Wxe^'^' + ^sign(a;)ue^l''l 
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and then, by using (3.36), we have 

1/(^)1' <II/IIhi(«)< / [2ul + {l + a)u']e-\y\dy<2C"'-. 
Jm 

Now we study the behaviour at infinity of the multipeakon solutions of the 
Camassa-Holm equation. 

Lemma 3.6. ('A-priori bounds) Let u be a multi-peakon solution to (1-1), with 
initial data u which belongs to the space defined in (3.35). Then for every 
f e K there exists a continuous function C{t), which depends on C"'" and on 
the energy E", such that 

• [[u^{t,x)+ulit,x)]e"^''Ux<C{t), (3.39) 
Js. 

• sup |PJ'(t,a;)|e"l^l < C{t) , (3.40) 

• \K\\LHn)<C{t). (3.41) 

Proof. Since \P^\ = P", it is sufficient to prove the second inequaUty with 
replaced bu P„. Setting 



I{t)= I [u^{t,x) + ul{t,x)]e''^''Ux, 



we want to achieve a differential inequality of the form 

p{t)<A + B.I{t), 

for some constants A and B which depend on the initial data u. We start the 
discussion proving a preliminary estimate for the function P". By applying the 
Fubini theorem to the identity 



dy (3.42) 



/ P''{t,x)e"\''Ux=^ [ e"l^lda; / 
Js. ^ Jr Js. 

we have to compute the following integral 

f e"kle-k-yl dx = ^^e-l^l + ^-^e"l«l for every yeR. (3.43) 
Jr 

For future use, we observe that the equality (3.43) holds for a G (—1, 1). Sub- 
stituting (3.43) in (3.42) and using the definition of the energy E'" we have 

/ P-{t, x)e'^\^\dx < — ^ + m . 

Js. (1 - a^) ^ ' 



3.2 A priori bounds 
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Having in mind the previous inequality, wc are able to estimate the time deriva- 
tive of the function /. From the equations (1.1) and (1.5) we have 



dt 



lit) = [ [2uut+{ul)^] e"l^lda; 



-2u{u^ + P^} + -(??), - (uul)^ - 2m^P« 

< -2 / (wP" + uul)^ e"l^l dx < -2u{P" + ul) e° 
Jr 

+2a / \u\iP'' + ul) e"\''\ dx 
Jr 



e"l^l dx 

oo 
-oo 



< 



2a|l?ilko 



2VE^ 



(£;« + 27) < Y^-j + 2/(i)] 



1-q2 

the previous inequality gives then a bound on the function /, that is 

I{t) < (C«'" + £"/2) exp i) • 

To achieve the estimate (3.40), set 

K{t)= P''(t,-)e"l-l 

Proceeding as before, fixed x GRwe compute the derivative w.r.t the time t of 
the function e~l^l * 



d_ 
dt 



-\x\ , '^x 

e ' ' * — 



I'm 
1 



2 ' ' 



-\--M{t,y)dy 
-v\ \ 



< ||u|UooP" + 



3 2 
u\\l 



uP" +mP" 



< ||u||i,ooP"(t,a;) + 



/ e-\-y\p^{t,y)dy 
Jr 

m I 



«kL-k-!/| 



dy 



< \\u\\L^P^{t,x) + ^e-«l^l||u||L=c /<:(<) 

in the same way, the derivative of e^'^' * is 

i (e-l-l * ^ ) < / e-l^-^l« (|P;| + \uu,\) dy 

2P«(t,a;)+ / e-\^-y\p''{t,y)dy 
\ Jr 



< ll^lliOO (2P"(t,X) 



1 — a- 



re-"l^'ii:(t) . 



Multiplying the previous two inequalities with e"l^l we get 



d_ 
It 



K{t) < 3 + 



^K{t) 
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which yields (3.40). 



To achieve the last inequality, we write 




< ( e-'^\y\dy+ ( ul{y)e°'\y\dy<-+I{t) 



Jr Jr 



where the last estimate is given by (3.39). 



□ 



3.3 Definition of the distance in the real hne 

In this section we define a metric in order to control the distance between two 
solutions of the equation (3) whenever we take initial data in the subspace Xa 
of H^{M.). It is constructed as in the spatially periodic case by resolving an 
appropriate optimal transportation problem. Let T = [0, 2ir] be the unit circle 
with the end points and tt identified. 




Figure 3.2: "Periodic" vs. "non-periodic" transportation map 

Consider the metric space (M^ x T,d^), with distance 

d'^{{x,u,Lo), {x',u',Lo')) = min{|a; — x'\ + \u — x'\ + \u> — 1} 
and for every function u G X^, let us define the Radon measure on x T 



for every Borel set A of x T. 

Deflnition 3.1. The set of transportation plans consists of the functions 
with the following properties: 




[l + ul{x)]dx 



3.3 Definition of the distance in the real line 
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1. ip is absolutely continuous, is increasing with its inverse; 

2. sup|a;-^(a;)|e"/2|^l < oo; 

3. /k |1 — i^'ix)] dx < oo. 




X 



Figure 3.3: Transportation plan. 

The conditions 2 and 3 are not restrictive. Indeed, thanks to the exponential 
decay of functions u,v G Xa, the measures cr" and ct" located on the graph of 
u and V respectively, have small mass at the infinity, and then a transportation 
plan which transports mass from one to the other can be almost the identity 
tp{x) « X (see fig. 3.3). In order to define a distance in the space X^, we 
consider an optimization problem over all possible transportation plans. Given 
two functions u, v in X^, we introduce two further measurable functions, related 
to a transportation plan ip: 

Mx)=snp{9e[0,l] s.t. 9-il + ulix))<{l + vl{^{x)))i^'ix)},i3.U) 
</,2(V'(a;)) =sup{6'G [0,1] s.t. 1 + uUx) < 9 ■ {l + vl{ip{x))) ip'{x)}. (3.45) 

The functions , (1)2 can be seen as weights that take into account the difference 
of the masses of the measure cr" and cr". In fact, from the definitions (3.44)- 
(3.45) one has 

(?ii(a;)(l + ul{x)) = <j)2{i^{x)){l + vl{il){x)))'ilj'{x) for a.e. x €R. 

According to the definitions, the identity max{^i(a;), (l)2{x)} = 1 holds. Altough 
the two measures (^icr" and (^2cr" have not finite mass, they satisfy (^icr"(^) = 
(l)2<7^{A) for every bounded Borel set A C x T. Thus, the functions and 4>2 
represent the percentage of mass actually transported from one measure to the 
other. A distance between the two functions u, v in Xa can be characterized in 
the following way. 

For every transport map tp G J^, let define X" = (a;, u(a;), 2 arctan'Ua;(a;)) 
and X" = {ip{x),v{il){x)),2 a.Tcta.nVx{tp{x))) and consider the functional 

J^{u,v)= [ rfO(X",X'')</.i(a;)(l+w^(x))da;+ / \l + uHx) - (1 + vli'iP{x)))iP' {x)\ dx . 

Jul Jul 

Since the above functional is well defined for every ip G J^, we can define 

J{u,v) = inf J'^{u,v). 
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The functional J here defined is thus a metric on the space Xa (see the previous 
Section 3.1). 



3.3.1 Comparison with other topologies 
Lemma 3.7. For every u, v G one has 

^- \\u-v\\li(m.) < J{u,v) < C ■ u||hi(k)- (3-46) 

Let (un) be a Cauchy sequence for the distance J such that C"''"" < Co for 
every n e N. Then 

i) There exists a limit function u e X^ such that Un ^ u in L°° and the 
sequence of derivatives Unx converges to in L^'(M) for p G [1,2[. 

a) Let be the absolutely continuous measure having density with re- 
spect to Lebesgue measure. Then there exists a measure ji whose absolutely 
continuous part has density such that /!„ H- 

Proof. The first inequality of (3.46) can be achieved by estimating the area 
between the two functions u and v. For every ip £ we can write 

\u — v\dx= \u — v\dx+ \u — v\dx 
Js. Jsi Js2 

where the two subsets Si and ^2 are 

• Si = {x : \x — ^{x)\ < 1 } = Uj[x2j-i,X2j], where in this union we have 
to take into account that these intervals may be either finite or infinite, 
possibly having Xj = ±oo for some j, 

• S2 = {x: \x-tlj{x)\ > 1}. 

The integral over 52 can be estimate in the following way: 

/ \u{x)-v{x)\dx<{\\u\\Loo + \\v\\Lo^) [ \x-^l;{x)\dx<{E"+E'')J{u.v). 

JS2 Jk. 

(3.47) 

The last inequality is given by the definition of the functional J. 

As far as the integral over Si is concerned, the intc;gral over 5*1 can be 
viewed as a sum of the area of the regions Aj in the plane M^, bounded by 
the graph of the curves u, v and by the segments with slope ±1 that join the 
points Qu{x2j-i) = {X2j-i,u{x2j-i)) and Qv{x2j) = {ilj{x2j),v{ip{x2j))), where 
{xi} = dSi. We have 



u{x) — v{x) \ dx < 2_\ ineas(Aj) . 



3.3 Definition of the distance in the real line 
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r ^ 








V 












Figure 3.4: L^— distance between two functions. 



The measure of the subset Aj is the area sweeped by the segment Qu{x) Qv{x). 
Recalling that in every set Aj the function ^ satisfies |a;| — 1 < \il){x)\ < |a;| + l, 
a bound on this area is given by 

meas(Aj) < / {\x - ii){x)\ + \u{x) - v{il))\)[l + ul + {I + vl{ii)))'4)']dx 
and then 

/ \u{x)-v{x)\dx < f i\x^i;ix)\ + \u{x)-v{i;)\)[l+ul + {l+vl{i;))il;']dx 

JSi JSi 

< J'^{u,v) + J''''\u,v) 

this inequality, together with (3.47), yields to 

||w-t;||i:,i(M) < C{u,v) J{u,v). 

Concerning the second part of the lemma, let us observe that even if the em- 
bedding of H^(M.) in L^(K) is not compact, the uniform exponential decay of 
the function u„ ((7"'"" < Co uniformly in n) allows us to extract a subsequence 
which converges to a function u in L^-norm and, by the unformly Hoder conti- 
nuity od such functions, \\un — u\\l'^ — > 0. This allow us to prove the property 
ii) by dealing with analogOTis arguments to tlios(^ developed in the periodic case. 
The proof of the second part of the lemma is thus perfectly similar to the one 
of the periodic case, once we take into account the exponential decay of the 
sequence u„. □ 
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3.3.2 Continuity of solutions w.r.t initial data 

Let uq and vq be two multipeakon initial data. The technique developed in 
Section 2.1 ensures the existence of two multipeakon solutions u{t),v{t) for 
(1.1) which conserve the energy unless interaction of peakon occurs. Suppose 
then that within a given interval [0, T] no interaction occurs neither for u{t) nor 
for v{t). The aim of this section is to prove the continuity of the functional J 
w.r.t. the initial data, namely we prove that there exists a continuous, positive 
function C{t) such that for t G [0,T] one has 

J{u{t),v{t)) < C{t)J{uo,vo). 

Lemma 3.8. Ifu(t) and v(t) are two multipeakon solutions defined in the inter- 
val [0, T] in which no interaction occurs, then there exists a positive, continuous 
function c{t) which depends only on the energies E^, E'" of the two solutions, 
such that 

^J{u{t),v{t)) < c{t)J{u{t),v{t)) for all t e [0, T]. (3.48) 

Proof. We compute the time derivative of the function J^{u{t),v{t)) with a 
particular choice of the transportation plan ij: = tp{t)- Given any G at 
every time t G [0, T] we construct by transporting the function ipo along 
the characteristic curves. More precisely, since no interaction between peakon 
occurs in the interval [0, T], the functions u{t, •), v{t, •) are Lipschitz continuous, 
then the flows ipl^, </3* solutions of the Cauchy problems 

^^iix) = u{t, ^i{x)) ip'iix) = X, 

^^Uy) = v{t,v=liy)) vl{y) = y, 

which are the characteristics curves associated to the equation (1.1), are well 
defined. Now, let a; e M. Vct) is defined as the composition 

V'(t)(a:) oV^o°K)"'(a:), (3.49) 

that is 

V'(t)(¥'t(y)) = ^*(^o(y))- 

The function tp(^f) belongs to and hence J'^'*) is well defined, in fact 

1. By the Property 1 of Definition 3.1 for the function -^o, and uniqueness of 
solution of ODE, the function is an increasing function. 

2. Let x € K. and (flj^{y) be the characteristic curve passing through x at 
time t. Evaluating \x — '(/;(t)(a;)|e"/^l'^l along this characteristic curve, and 
computing the derivative w.r.t. t we obtain 

||^t(y)-^*(^o(y))|e"/^l^"(^)l 



< 



\u{t,x)-v{t,i)^t){x))\ + ^\u{t,x)\ ■ \x - llj(^t){x)\ 



^a/2\x\ 



3.3 Definition of the distance in the real line 



57 



by properties (3.38), (3.39), and since u, v are Lipschitz continuous in 
[0,T], there exists two functions ci{t), C2{t) such that 

||a;-^(,)(a:)|e«/2H < ci(i)k - V(t)(a;)|e«/2|-I + C2(i) 

by Gronwall Lemma and the hypothesis \x — ?/;o(.T)|e"/^l^l < Co, the pre- 
vious inequality gives the Property 2 of Definition 3.1 for tp(^f^ 

|x-V.(t)(x)|e"/2W <Ci(t) = (co + J\2{s)ds^ e/o-iW<i^ (3.50) 



3. The last property can be achieved by choosing the change of integration 
variable x = ip*^ (y) 

/ |l-V(t)(a;)|da; = / \1 - ij^,^{x)\{v>iy {y) dy 

= J \{'fiiny)-{'pl)'{MyWo{y)\dy 

< I My{y)-l\dy+ [ \{^iy{y)-l\dy 

+ / \l-i;'o{y)\dy. 
Jr 

Since 

My{y)-1\< f\u,{s,x)\-\{<piy{y)-l\ds+ f\u,{s,x)\ds 
Jo Jo 

(and a similar estimate for (^*) and Ux,Vx & L°° , by the Gronwall lemma 
the first two integrals of the previous formula are bounded by an absolutely 
continuous function C{t) in the interval [0, T] and then also Property 3 of 
Definition 3.1 holds. 

At the transportation plan V'(t) we associate the functions <i)i \ 4>2^ defined 
according to (3.44), (3.45), the functional J'^M is thus 



J^(*^{u{t),v{t))= / d'>{X''{t),X^t))4'\x){l+ul{x))di 



I l + "x(a^)-(l + ^^x(V'(t)(a^)))V'(t)(a^) 

Jr 



dx. 



By deriving J'^w {u{t), v{t)) w.r.t. t and computing the change of variables along 
the characteristics, the previous derivative can be estimate by the sum of the 
following terms (we leave out the dependence on the integrable variable when 
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it is not essential) 

h = [ \u{t,x)-v{t,^^t){x))\(t>'^\x){l+ul{t,x))dx 
Jr 

< (1 + \\u{t)\\Loo + \\vit)\\Loo)J^^*->{u{t),v{t)) , 

72= / \P:{t,x)-P^{t,i,^t){x))\ct>^^\x){l + ul{t,x))dx, 
Jr 

2u\t)-ul(t)-2P-{t) 2v^{t,i;(t))-vUt,i^(t))-2P''{t,i^(t)) 



h 



Jr 



l + ul{t) 
b^^\l + ul{t))dx, 



the term due to the variation of the base measure 

74 = 2 / d^(X"(t),X''(i))-M,(<)(u2(t)-F"(t))dx. 
Jr 

and the terms due to the variation of the excess mass 



dx . 



Let us start to estimate the term I2. By definition, the difference of P" and 
is written in convolution form 



/ \e 1^ ^1 sign(x - y) 
Jr I 



dy 



^(t){y) 



by this inequality, we can estimate the term I2 by the sum of the following 
integrals 



A = [ il + ul{t,x)) [ e-\--y\\u\t,y)-v\t,y)\dydx 
Jr Jr 



B 



{l+ul{t,x)) 



C = I {l + ul{t,x)) f 
Jr Jr 



e 1^ y\s\gn{x-y)[v'^{t,y)-v'^{t,ij(t){y)W(t){y)]dy 

V {t,^(t){y)) + — 



dx 



D 



1 



-\x-v\ 



sign(x -y)- e-l^w(-)-^(*)(^)l sign(V(t) (a;) - • V(t)(2/) dydx 



{l + ul{t,x)) / e-^--y\sign{x-y)[ul{t,y)-vl{t,ilj^t^{y))i;[t){y)]dy 

Jr 

A. Switching the order of the two integrals, the term A is bounded by the 

L^-norm of the difference between u and v: 

A <{\\u\\loo + \\v\\loo) [ \u{t,y)-v{t,y)\ [ e-^^'-y^il + ul{t,x)) dxdy 
Jr Jr 
< (2 + 7;")(||u||loo + ||i;||Loo)||u(i) - v{t)\\Li 



dx 
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and then, by Lemma 3.7, A < C{u,v)J{u{t),v{t)). 
B. Define 



F{y)= f {v^{z)-v\i^^^,){z))i,[,^{z))dz= r v\z)dz 



we have, integrating by parts 



/ 



sign(a; - y)F'{y) dy 



;2|F(a;)|+ / e-l--^l|F(y)| 
Jr 

< Mloo\x - i,{x)\ + [ e-\''-y\\y-^{y)\dy 
Jr 

moreover, substituting the previous expression into the term B we obtain 

B<2E^J^^^^{u{t),v{t))+ [ {l + ulit,x)) f e-l--^l|y-V(t)(y)|dyrfa; 

Jm JR 

= 2E''{j^iHuit),vit))+ [ |2/-V(t)(2/)| [{l+ul{t,x))e-^^-yUxdy] 
I Jr Jr 

< 2E^{3 + E^) ■ J'l'^*\u{t), v{t)). 



C. Observe that since the fimction y '<l\i){y) is non decreasing, the quantities 
X — y and ^(j)(a;) — '4'{t){y) have the same sign, and since the function t ^ e~l*' 
is Lipschitz continuous either in (— oo,0) or in (0, +oo) we have 



-W-v\ _ p-\i'{t){x)-iint){v)\ 



< g-min{|.-y|,|V'(*,W-^w(y)|} \\x - y\ - | ^'(t) (x) - (?/) | | 

< g-min{|x-yMV.(,)(x)-V>w(2/)|} _ ^^^^ | + |y _ ^^^^^^^Q 



now 



-min{|a;-y|,|V'(t)(a;) - < -\x - y\+2Ci{t), 

where C\{t) is the function (3.50), related to the Property 2 of Definition 3.1 
V'(t), then 



C< e^^W/ \y-^{y)\ 



v\t,i,(t){y)) 



■ [ {l + ul{t,x))e-^''-yUxdy + +2E^ [ {1 + ul{t,x))\x - tlj^t){x)\dx 
Jr Jr 

< \{2 + E'')e'^'^*\l + \\u\\Loo) + 2E^ J^w 

D. Here we can use the estimate given by the change in base measure. Since 

/ l + ul{t,x)-{l+vl{t,i^(t){x)))^P[t^{x) dx<J^^^>{u{t),v{t)), 
Jr 
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we obtain 

D<1 [ {l + ulit,x)) [ e-l^-^l il + ul{t,y))-{l + vl{t,^l;^t)iy)W[,^{y) 
{l+ul{t,x)) 



dy dx 



e 1=" y^sign{x-y) V'(t)(y) - 1 



< 



+ ul{t,x)) 



(V'(t)(a;)-a;) 



-\x-y\ 



+ (l + £;")J'^w {u{t),v{t)) 
< 2{2 + E^)J'l'w{u{t),v{t)) 



/ {^{t){y)-y) 



dy 



dx 
dx 



where in the last estimate we integrated by part as in the term B. 

The control for the terms /s, /4 and can be obtained exactly as the ones in 
[10], whom we refer the reader to. The previous estimates implies that there 
exists a smooth function C = C"''" {t) which depends only to the variable t and 
to the initial data u, v such that 

f^J^{u,v)<C'^^^t)J^{u,v) 

which yields 

J{u{t),v{t)) < J(u(s),u(s))el-''»*'^''''('')H for every s,i e R. 



□ 



3.4 Proof of the main theorems 

Thanks to the analysis in the previous sections (3.1 for the spatially periodic 
case, 3.2 and 3.3 in the whole real line), we now all the ingredients toward a 
proof of Theorems 1.1 and 1.2. Since the two cases are analogous, we prove them 
in the periodic case. The estimates in (1.14) follow from Lemma 3.1. Given an 
initial data u e H^^j., to construct the solution of the Camassa-Holm equation 
we consider a sequence of multi-peakons Un, converging to u in -ffper- Then we 
consider the corresponding solutions t Un{t), defined for all n > 1 and t gM.. 
This is possible because of Lemmas 2.2 and 2.3. 

We claim that the sequence M„(i) is Cauchy in L^^j. . Indeed, by Lemma 3.1 
and Lemma 3.3, 

\\Umit) - Un{t)\\ I < C ■ J(um{t), Un{t)) 
-'-'per 

< C ■ e-\'\ J{uUO), «n(0)) < ■ e«l*l||u„(t) - u„{t)\\^ 

Therefore, Un(t) u(t) in L^^^, for some function u : M -^pcr • By interpo- 
lation, the convergence Un ^ u also holds in all spaces ipgj., 1 < p < oo. The 
continuity estimates (1.15)-(1.16) now follow by passing to the limit in Lemma 
3.3 and 3.4. 
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It remains to show that the hmit function u(-) is actuaUy a sohition to the 
Camassa-Holm equation and its energy E{t) in (1.13) is a.e. constant. Toward 
these goals, we observe that all solutions u„ are Lipschitz continuous with the 
same Lipschitz constant, as maps from K into L^^j. . Indeed 

lkr»,t|k|„ < ||Mn||L~ • ||Mn,cc||z,|„ + 

As a consequence, the map t u{t) has uniformly bounded iJp^r norm, and is 
Lipschitz continuous with values in L'^^^. . In particular, u is uniformly Holder 
continuous as a fimction of t,x and the convergence Um(t,x) u{f.,x) holds 
uniformly for t in bounded sets. Moreover, since L^^^. is a reflexive space, the 
time derivative Ut{t) G L^^^. is well defined for a.e. t gM.. 

We now observe that, for each n > 1, both sides of the equality 

^Un = -Un Un,x " -P"" (3.51) 

arc continuous as functions from M into L^^^ , and the identity holds at every 
time t G K, with the exception of the isolated times where a pcakon interaction 
occurs. 

At any time t where no peakon interaction occur in the solution u„, we 
define /Zj"' to be the measure with density w^(i, •) + \un,x{t, ■) w.r.t. Lebesgue 

measure. By Lemmas 3.3 and 3.2, the map t ^ fi^^^ can be extended by weak 
continuity to all times t e M. We can now redefine 

(3.52) 

where fit is the weak limit of the measures Because of the convergence 

J(un{t), M(t)) — > 0, by Lemma 3.2 the map i > /it is well defined and continuous 
w.r.t. the weak topology of measures. Using again Lemma 3.2, we can take the 
limit of (3.51) as n — > oo, and obtain the identity (1.12), for every t G R and 
P = P« defined by (3.52). 

For each n, the total energy /Xj"^(]0, 1]) = E^^ is constant in time and 
converges to E"^ as n — > oo. Therefore we also have 

Mt(]0,l]) =£" = / [u'^{x)+ul{x)]dx ieR. 
Jo 

To complete the proof of Theorem 1.1, it now only remains to prove that 
the measure is absolutely continuous with density 

w.r.t. Lebesgue measure, for a.e. time t gM.. 



1 



L2 
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In this direction, wc recall that, by (1.5), the function w = u^^ satisfies the 
linear transport equation with source 

Moreover, along any characteristic curve 1 1-^ ^(t) by (1.3) one has 



di 



2 arctanu„_i: (t, ^{t, x)) 



1 + Ul^x 



< - 2 ' (3-53) 



whenever ^ is sufficiently large. For £ > small, consider the piecewise affine, 
27r-periodic function (see fig. 3.5) 



f{0) = { 



6 


if 


< 6* < 1 


1 


if 


1 < 61 < TT - £, 




if 


TT — £<0<7r + e. 


-1 


if 


7r + e<6l<27r-l, 


6l-27r 


if 


2tt-1<6<2'k. 




and define 



'n{t) = / ip{2BxctanUn,x{t,x))u'^^^{t,x)dx . 
Jo 



By (1.3) and (3.53) we now have 

tA^^^^tJ ul,{t,x)dx-C- [ 

^£ J{2arctanu„ j,e[7r-£,7r]U[-7r,-7r+£]} J{2 arctanM„ ;,e[-l,l]} 

(3.54) 

for some constant C, independent of s,n. Since all functions /3„ remain uni- 
formly bounded, by (3.54) for any time interval [r, r'] we obtain 

/ / ul,^{t,x)dxdt<eC' ■ {l + T -t), (3.55) 

J T J {2 arctanun,x£['n"— £,7r]U[— TT,— TT+e]} 
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where the constant C depends only on the H^^,^ norm of the functions ?i„, hence 
is uniformly valid for all n, e. Because of (3.55), the sequence of functions ^ 
is equi-integrable on any domain of the form [r, r'] x [0, 1]. Namely 

lim / / ul^^{t,x)dxdt = 0, (3.56) 

i^^°°Jt ^{xe[o,i], < ,^>k} ' 

uniformly w.r.t. n. By Lemma 3.2 we already know that 1 1 ^ (t) — (i) 1 1 i — > 

for every fixed time t and 1 < p < 2. Thanks to the equi-integrability condition 
(3.56) we now have 

ul^^ul in L1([t,t']x[0,1]). 

By Fubini's theorem, this implies 



lim / ui {t,x)dx = I ui{t,x)dx 



for a.e. t G [t, r']. At every such time t, the measure /Zf is absolutely continu- 
ous and the definition (3.52) coincides with (1.2). This completes the proof of 
Theorem 1.1. □ 



3.5 Uniqueness 

Before proving Theorem 1.3, we remark that the solution satisfying all condi- 
tions in Theorem 1.1 need not be unique. 



Example 3.1. Let u = u{t,x) be a solution containing exactly two peakons 

of opposite strengths pi{t) — —p2{t), located at points qi{t) = —q2{t) (sec fig. 
3.6). We assume that initially pi{0) > Ogi(O) < 0. At a finite time T > 0, the 
two peakons interact at the origin. In particular, as f — > T~ there holds 

pi{t) ^ 00, P2{t) ^ -oo, qi{t)^0, q2{t)^0. 

Moreover, ||it(t)||^^ 0, while the measure fit approaches a Dirac mass at the 
origin. We now have various ways to extend the solution beyond time T: 

u{t,x) = , 

u{t, x) = -u{t - T, -x) , (3.57) 

Clearly, u dissipates all the energy, and does not satisfy the identity (1.13). 
The function u in (3.57) is the one constructed by our algorithm in Section 2.1. 
However, there are infinitely many other solutions that still satisfy (1.13), for 
example 

■u(t, x) = u(t, X — b) 

where u is as in (3.57) and 6 7^ 0. The additional condition in Theorem 1.3 rules 
out all of them, because as r ^ T+, the corresponding measures jlr approach 
a Dirac mass at the point x = b, not at the origin. 
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Figure 3.6: Two solutions of the peakon-antipeakon interaction 



We can now give a proof of Theorem 1.3. As a first step, we extend our 
distance J to a larger domain V, consisting of couples {u,^), where u € -ffper 
and /X is a positive (spatially periodic) measure whose absolutely continuous 
part has density + w.r.t. Lebesgue measure. This extension is achieved 
by continuity: 

J({u,ii), = liminf J(u„,u„) 

where the infimum is taken over all couple of sequences (w„,u„)„>i such that 
\\un- u\\l<=° ^ , - wIIloo ^ , 

2 ~2 ~ 

We observe that the flow $ constructed in Theorem 1.2 can be continuously ex- 
tended to a locally Lipschitz continuous group of transformations on the domain 
V. 

Now let 1 1— > u{t) be a solution of the Cauchy problem (1.1), (1-11), satisfying 
all the required conditions. In particular, the map t ^ [u{t),jlt) is Lipschitz 
continuous w.r.t. the distance J, with values in the domain V. 

Calling t {u{t), fit) = $*(«, u^) the unique solution of the Cauchy problem 
obtained as limit of multi-pcakon approximations, wc need to show that u{t) = 
u{t) for all t. To fix the ideas, let t > 0. By the Lipschitz continuity of the flow, 
we can use the error estimate 

j{(u{t),fit) , (M(t),Mt)) < e^^* J limmf~j(^{u{T+h),jlr+h) , ^h{u{r), f'r)) dr 

(3.58) 
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For a proof of (3.58), sec [6, pp. 25- 27]. The conditions stated in Theorem 1.1 
now imply that, at almost every time r, the measure jit is absolutely continuous 
and the integrand in (3.58) vanishes. Therefore u{t) = u{t) for all t. □ 

We can now prove that, in multi-peakon solutions, interactions involving 
exactly two peakons are the only possible ones. 

Corollary 3.1. Let t i— > u{t,-) he a multi-peakon solution of the form (1.6), 

which remains regular on the open interval ]0,T[ . Assume that at time T > 
an interaction occurs, say among the first k peakons, so that 

^lim qi[t) = q i = l,...,k. 

Then k = 2. 

Proof. We first observe that the Camassa-Holm equations (1.1) are time re- 
versible. In particular, our proof of Theorem 1.3 shows that the solution to a 
Cauchy problem is unique both forward and backward in time. 

Now consider the data {u{T), ^t) & 'D, where /xt is the weak limit of the 
measures fit having density ?i^(t) + u^(t) w.r.t. Lebcsguc measure, as t ^ T— . 
By the analysis in Section 2.1, we can construct a backward solution of this 
Cauchy problem in terms of exactly two incoming peakons. By uniqueness, this 
must coincide with the given solution u(-) for all t G [0,T]. □ 
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Part II 



The discrete Boltzmann 
equation 



Chapter 4 



Symmetry groups of 
differential equations 

In this chapter we introduce the theory of the symmetry groups apphed to 
differential equation, which is a tool that will fits in the study of evolutionary 
equations. The goal of this Chapter is to develop a useful method that will 
explicitly determine the symmetry group for the system of discrete Boltzmann 
equation, which will be the starting point of the discussion of Chapter 5 for 
the blow-up issue. The key point is to transform the equation which has an 
asymptotic blow-up at a time T into an equation, related to a rescaling of the 
first equation, which approach a steady state as r goes to infinitive (see Section 
4.2). 

4.1 Group and differential equations 

The symmetry group of a system of difi^erential equations is the largest local 
group of transformations acting on the independent and dependent variables 
with the property that it transforms solution of the system to other solution. 
In the first part of this section we review a general computational method for 
(almost) any given system of differential equations. For more information about 
the application of group theory to the differential equation, we refer the reader 
to [42, 43]. 

We start recalling some useful definition in the abstract theory of Transfor- 
mation Group. 

Definition 4.1. Let M be a smooth manifold. A local group of transformations 
acting on M. given by a (local) Lie group G is the couple ('^,*) where 

• ^ is an open subset {l} xMC^CGxM. 

• a smooth map ^ -.^ ^ M. 
satisfy the properties 
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(a) If {h, x) G"^, {g, ^{h, x)) G and {g-h,x) G'^ then 

<l!{g,^{h,x)) = ^{g-h,x). 

(h) For all x G M, *(t, x) = x. 

(c) If{g,x)e'^ then *(.g, .t)) G anrf *(.g-\ a;)) = a;. 

For brevity, when it does not make confusion, we denote '^{g,x) by g ■ x. 

Definition 4.2. Let G be a local group of transformation acting on a manifold 
A4 . A subset S C ^A is called G— invariant, and G is called a symmetry group 
of S if whenever x G S and g G G are such that g ■ x is defined, then g ■ x G S. 

Remark 4.1. In our applications, as far as the difTcrcntial equation is con- 
cerned, the subset <S will be usually the graph of the solution of the differential 
equation 

^x,u,...,D"u,...) =0, (4.1) 

where i.e. set of solutions determined by the c;oninion zeros of collection of 
smooth functions $ = ($i, . . . , where $j = ^j{x, u, . . . . . . ) depends 

on the variables x and the unknowns and their derivatives u, D"u, . . . , and 
where, for every multi- index a = {ai, . . . ,am), indicates the differential 
operator 



dxi ) \ dx 

In this context, it is thus useful to introduce the graph of a function u : O — > i?" 
defined on a open set ^ C M"* 

Tu = {{x,u{x)) -.xgQ.} 

which is a smooth submanifold of R™ x E" . The action of a given transformation 
g G G maps the graph r„ into the subset g-Ty, = {{x, u) = g- {x, u) : {x, u) G Ty} 
which is not necessarily the graph of a function u. However, since G acts 
smoothly and the identity i G G leaves r„ unchanged, by restricting the domain 
f2 for every g G G near the identity the transformation £f • r„ is the graph Ty of 
a function u. 

As an example of action on the graph of a function, let consider a vector 
field v(^, ry) : M = — > which can be seen in local coordinate 

{x \ , . . . , Xmi ^1 ) • • • ? '^n) 

i=l j=l 

it acts on a smooth scalar function cj) : K™+" — > M as a derivation 

(t>{x + s^,u + er])-(l){x,u) 
V • (j)(x,uj = \im . (4.2) 

s—>o e 



4.1 Group and differential equations 



71 



The most important operation on vector fields is their Lie bracket or commuta- 
tor. Whenever we think that two vector field v, w act as a derivation, their Lie 
bracket [v, w] is the unique vector field satisfying 

[v, w] • = V • (w •(/)) — w • (v • (/)) for all smooth functions (j) (4.3) 

The integral curve of the vector field v is a smooth parametrized curve P{9) = 
(x, u) whose tangent vector at any point coincides with the value of v at the 
same point: 

l^ = ^(^)- 

Starting from a given initial data -P(O) = P = [x, u) the corresponding integral 
curve is often denoted by the suggestive exponential notation 

exp(6'v)P. 

From the existence and uniqueness of solution to systems of ordinary differential 
equations we easily obtain the semigroup property for the flow generated by v: 

exp(Ov)P = P (4.4) 
exp[(6ii + 6i2)v]P = exp(6iiv)[exp(6l2v)P]. (4.5) 

From these formulas, compared with the property (a)-(b) of definition 4.1, we 
see that the flow generated by a vector field is the same as a local action of 
the Lie group R on the manifold M. which is called a one-parameter group of 
transformations. The vector field v is called the infinitesimal generator of the 
action, v is also called an infinitesimal symmetry generator for (4.1) if the map 
P ^ exp(6'v)P transforms the graph of a solution u into the graph of another 
solution. 

Remark 4.2. Recall that if 9 is sufficiently small, there exists a neighborhood 
V such that the set 

Fe = {exp(0v)(.T,7i) : {x,u) e F„} 

coincides on V with the graph T^e of a smooth function . By using a Taylor 
expansion, the flow exp(0v) maps the point (x, u) into the point 

exp(^v)(x, u) = {x + 9^{x, u) + o{0), u + er}{x, u) + o{9)) 

therefore, the differentiation w.r.t. 6 at the origin yields the useful formula 



= r]- Wu{x) ■ ^. (4.6) 

61=0 



The previous formula gives a necessary condition in order to prove that a 
particular vector field v = {^,r]) is an infinitesimal symmetry generator for the 
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difFcrcntial equation (4.1). If the function another solution to this equation, 
then $(a;, u^, . . . , D"u^ , . . . ) =0. Differentiating w.r.t. in we get 



61=0 



(4.7) 



In the following we shall prove that the previous condition is sufficient in order 
to construct an infinitesimal symmetry generator v. 

Proposition 4.1. Let G he a connected group of transformation acting on the 
manifold M . A smooth real-valued function C, : M. ^ is an invariant function 
for G if and only if 

V • C = for allxeM (4.8) 
and every infinitesimal generator v of G. 

Proof. Suppose that C, is an invariant function for G. According to (4.2), if 
x&M 

— C(exp(6'v)x) = V • C[exp(6'v)x] 

since C, is invariant, setting = it proves the necessity of (4.8). Conversely, if 
(4.8) holds then ^(exp(0v)a;) is a constant for the connected subgroup {exp((?v)} 
oi Gx = {g & G : g ■ X \s defined}. But by the properties of the Lie group, every 
element of Gx can be written as a finite product g = exp(0-'^Vij) • • • exp((?'^Vi^ ) 
for some infinitesimal generator Vi of G, hence C,{g ■ x) = C,(x) for all g € Gx- 

□ 

In a similar way we can prove the following theorem which gives an infinites- 
imal criterion of invariance for a general equation 

<I>(x) = X e M 

that will be useful whenever we are concerning a differential equation 

^{x,u,...,D°'u,...) =0. 

Theorem 4.1. Let G be a connected local Lie group of transformations acting 
on a p— dimensional manifold M . Let $ : — !■ M' , I < p, define a system of 
equations 

$,(a;) = 1^=1,..., I (4.9) 

and assume that the system has maximal rank at every solution x of the system, 
namely 



rank 



^(x) = 0. 



Then G is a symmetry group of the system if and only if 

v.$^(a;)=0 v = l,...,l, $(a;)=0 
for every infinitesimal generator \ of G. 



(4.10) 
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Proof. Let xq be a solution of the system (4.9). As in Proposition 4.1, the 
necessary condition follows from differentiating w.r.t. e the identities 

$^(exp(£v)xo) = 

and setting e = 0. Conversely, by using the maximal rank condition, we can 
choose local coordinates y = {y^, . . . , y™) such that xq ~ and $ has the simple 
form = {yi, . . . ,2/')- ^ be an infinitesimal generator of G, which can 
be expressed in the new coordinates as 

The condition (4.10) turns to be 

viy") = C{y) =0 for all = 1 ... Z 

whenever y^ = ■ ■ ■ = y^ = 0. Since the flow 0(^) = exp(^v)a;o satisfies the 
system of ODE 

^;.*(o) = o i-i...m, 

the uniqueness of the solution yields to conclude that ([>' (0) = ior 6 sufficiently 
small. exp(^v)xo is thus again a solution to ^{x) = 0. As in Proposition 4.1, 
by the properties of the connected local Lie group G we gain the result. 

□ 

The previous theorem can be adapted for the differential equation in order 
to get sufficient condition for obtain an infinitesimal symmetry generator. The 
equation (4.1) contains not only the unknowns u but also its derivatives, so 
in order to use Theorem 4.1 wc can think that the solution is a point which 
contains all of these functions. To do this we need to prolong the basic space 
representing the independent and dependent variables under consideration to 
a space which also represents the various partial derivatives occurring in the 
system. 

If we consider function u : O C — > M", the number of derivatives of order 
k is 

n + k-1 
k 



Pk 



If N is the maximum order of the derivatives involved in the differential equation 
(4.1), we introduce thus the N-thjet space ilxU^^'' = nxM"' xW^ x ■ ■ -xW, 
whose coordinates represent all the derivatives of the function u from to N. 
If u is a function whose graph lies in a manifold A4 C CI x M™, we define its 
prolongation 

prW« = {u, (D"i«)|,,|=p,, . . . , (i5""w)|,„|=p„) 
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whose graph hcs in the A^— th jet space M^'^^ = M x x • • • x RP". From 
this point of view, a smooth solution of the given system (4.1) is a function u{x) 
such that 

prWM(x)) = ;/ = !,...,/, 

whenever x Hes in the domain of u. It means that the graph of the prolongation 
of u must lie entirely within the subvariety of the zeroes of the system. 

Now suppose that G is a local group of transformations acting on A4. The 
prolongation of G as a local action group on the prolonged manifold Ai^^^ 
is defined so that if g G G, transforms the derivatives of u into the 

corresponding derivatives of the transformed g ■ {x,u). To evaluate the action 
of the prolonged pr^^^g on a couple {xo,u^^) we simply choose a particular 
function / whose derivatives agree, up to A^— th order, to the point {xo,Uq^^), 
apply the action g to / and then prolong g ■ f. Last, we have to define also the 
prolongation of a vector field up to the order N . It follows by viewing it as the 
infinitesimal generator of the corresponding action group pr(^)[exp(0v)]: 



(N) ■ <^ 



pr(^)[exp(^v)]. 



Writing 

A/ ^ 9 , , (9 

V = ^{x,u)—+g{x,u) — , 

the general formula of such a vector field is given by [42, Theorem 2.36] and it 
is the following formal expression 

prWv = v+ ^ (4.11) 

l<|a|<Ar 

where 

(f>" = D'' (r/-^V^^)+^ Vu". 
For future use, the prolongation of the Lie bracket vector field is 

prW[v,w] [pr^v, pr^w]. (4.12) 

By applying Theorem 4.1 to the equation ^{x, pr'^^w(a;)) = we obtain the 
following theorem, which agree with the formula (4.7) 

Theorem 4.2. Suppose 

$^(x, prWu(a;)) = u = l,...,l, 

is a system of differential equations of maximal rank defined over M. C a; x K". 
If G is a local group of transformations acting on A4, and 



pr 



(JV) 



V • ^„{x, pr^^^u) = 0, for all v = I, . . . ,1, 



for every infinitesimal generator v of G, then G is a symmetry group of the 
system. □ 
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4.2 Symmetries and blow-up 

In this section we enter in deep detail of the application of local symmetry 
groups to the blow-up issue. For a more complete description, the reader can 
see [5] . Our goal is to apply the theory previously developed in order to choose 
an appropriate infinitesimal symmetry generator which describes the asymptotic 
behaviour of a solution which has blow-up in finite time. 

We shall consider an evolution problem on a Banach space {\\ ■ \\,E) 

X = fix). (4.13) 

Assume that there exists trajectories that blow-up in finite time 

lim \\x(t)\\ = +00. 

Suppose that there exists a second vector field g such that 

• Trajectories oi y = f{y) — g{y) do not blow-up. Instead, they approach a 
steady state y as time goes to infinity. 

• There is an explicit computable transformation that maps a trajectory 

s I— > y{s), s G [sq, +oo[ into a trajectory 1 1-^ x{t), t € [to, T[. 

In this case we could first accurately study the asymptotic behaviour of y{s) 
as s ^ -l-oo, and then recover information on the behaviour of the blowing-up 
solution x{t). To implement this approach, it is clear that the auxiliary vector 
field g must be carefully selected. 

According to the theory developed in the previous section, we shall look for 
condition on 

^x«=5(x«) (4.14) 

in order to have that = exp{9g)u gives another solution to the equation 
(4.13), provided that we have existence and uniqueness of solution to the Cauchy 
problem (4.14). 

As the following analysis will show, the crucial assumption on g is the relation 

f + [f,g] = 0. (4.15) 
Lemma 4.1. Let g be a vector field such that 4-15 holds. Then the vector field 

is an infinitesimal generator of symmetry group. In other words, if t^ x{t) is 
a solution of (4-13) then 

{exp(0v)(t,a;(t)) : t G /} 
is the graph of another solution 

x\t) =exp{9g){x{e''t)). 
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Proof. To check that is indeed a solution, we write 

^=e' Jac [expieg]] f (exp(-eff)a;^) . 

We claim that the right hand side of the previous formula coincides with /(x^), 
and this can be done by proving 

v{e) = Jac [eMOg)] f (exp(-^ff)t/) = e-^f{y) for all y G E. (4.17) 

Trivially, i;(0) = ,f{y). As far as the derivative of v is concerned, let use the Lie 
bracket property (see [25]) 

[/ g] = lim [exp(£g)] g(exp(£/)) - g 
s—>o e 

and the hypothesis (4.15) 

dv{0) ^ Jac[exp[(£+e)]g)]/(exp[(-e-£)g]i/)- Jac[exp(gg)]/(exp(-eg)i/) 

= - Jac [exp(6ifif)] • [g, f] exp{-eg)y = -v{e). 
The function v{6) is thus 

v{e) = e-%(0) = e-'fiy). 

□ 

Theorem 4.3. (Blow-up rescaling) Consider two vector fields f,g satisfying 

(4-15). Let y : [0, (X)[i-^- E be a solution to 

y{s) = fiy{s)) - g{y{s)). 

Then the function x : [0,l[i-^ E defined by 

x{t) = exp{sg){y{s)) s = In ^ (4.18) 

is a solution of (4-13). 
Proof. Notice that 



g{ex.p{sg)y) = ■^exp{eg){exp{sg)y) 



= Jac[exp(s5)]s'(y) 

e=0 



Let compute the derivative w.r.t. s of the function exp(s(/)t/(s). By the previous 
identity and by (4.17) we have 

^ exp(s5)y(s) = g{exp{sg)y{s)) + Jac[exp(sp)]2/(s) 

= 5(exp(sg)y(s)) + Jac[exp(s5)]/(?/(s)) - Jac[exp(s5)]s'(y(s)) 
= Jac[exp(ss()]/(exp(-s5)a;) 
= e-^f(a;) 

now we can check that (4.18) is a solution to (4.13). 

^x{t) = ^ exp{sg)y{s) ■ ^ = e-'f{x) ■ e' = f{x). 

□ 
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Figure 4.1: Connection between y{-) and x{-) 



We shall apply the above theory to a special case of partial differential equa- 
tions. Inside the equation (4.1) we highlight the time variable x {t,x), ob- 
taining the evolution equation 

F{x, pr^^^i) - = F{x, u,..., D"u, ...)-ut = (4.19) 

where the derivatives on involves only derivatives with respect to the spatial 
variable x. In the following, we indicate with [u] the prolonged function u in 
the space Ai x U'^^^ Thus, the total derivative of a function F[u{t)] takes the 
form VF ■ [ut], where [ut] = (w*, . . . , • • • )• 

The main tool is Theorem 4.2, which gives characterization on the symmetric 
vector field. Prom this theorem we obtain 

Theorem 4.4. An evolutionary vector field 

v = -t|+^(:,,.)^ + <^(.,.)|- 

is an infinitesimal symmetry generator for (4-19) if and only if for every function 
u one has 

{F+[F,G])[u]=0 (4.20) 

where 

F[u] = F{x, pr^^^w), G[u] = -^{x, u)Vu + (jj^x, u), 
[F,G][u] = hm G[n + eFM-G[u]-F[u + eGM-,F[u]_ 

£—>0 s 

Proof. Given any smooth solution m, let 6* sufficiently small and be the trans- 
formed function by the vector field v, according to Remark 4.2. As Theorem 
4.2 states, the thesis follows once we prove (4.7). By formula (4.6), substituting 
(^,r/) with ((— ^),(/>) we have 



du^ , ^ ^ , 

= + tut — z, ■ vu = tut + G[u\. 

0=0 



dO 
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Thus, formula (4.7) becomes 



VF[u] ■ [tut + G[u]] - F[u] - tVF[u] ■ [ut] - VG • [u]t 
-F[u\ + VF[u] • [G[u]] - VG[u] ■ [F[u]]. 



thus 

d^lu'^ 



= 



dO 

if and only if [F, G] [u] = -F. □ 

4.3 A group of symmetry for the discrete Boltz- 
mann equation 

The work plan presented in the previous two sections fits in the study of the 
blow-up rate of discrete Boltzmann equation. Let consider the system of PDE 
(see [41]) 

dtUi + Ci'VxUi = '^aijkUjUk i = l,...,l (4.22) 

jk 

where 

• {t,x) & Rx R^, 

• Cj e M*^ plays the role of velocity, 

• Ui{t,x) is the density of particles having speed Cj, 

• (lijk &re the coefiicients of the quadratic collision term, with ciijk = if 
j = k. 

We implement Theorem 4.4 in order to recover an evolutionary vector field which 
generates a group of symmetry for the system (4.22). Looking for a vector field 
of the form 

-*|+^(^)-V. + 0(.)-|^ 

where 



d d d \ d_ ^ /_d_ _d_ _d_ 

dx\ ' dx-2 ' dxz / ' du \ dui ' du^ ' ' ' ' ' dui 



we have to find condition for which (4.20) holds. Note that in this case the two 
maps F and G are 

{F[u])i = -Ci-Va;+^ aijkUjUk {G[u])i = -^{x)-'VxUi+(j)i{u) i = l...l. 
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Theorem 4.5. Setting 

^{x) = -X (j)i{u) = Ui, 

the identity 

[F,G][u]+F[u]=Q 
holds for every u solution to (4-22). 
Proof. Let compute (4.21). 

{G[u + eF[u\] - G[u])i _ x-\7a^(ui- sa ■ V^Ui + e J^j^ aijkUjU^ 

Ui - eCi ■ VxUi + S^jk O-ijkUjUk - X ■ VxUi - Ui 



+- 



= -x-Wx (ci • WxUi) + {F[u])i 

jk 

(4.23) 

{F[u + sG[u]] - F[u])i _ -Ci ■ VxUj - sCiVxjx ■ VxUi) - eciV^Ui 

e e 

J2ik (lijkiuj + ex ■ VxUj + euj){uk + ex ■ V^Uk + euk) 
H 

-CiVxUi + J2jk aijkUjUk 

e 

+ ^^aijk{2ujUk + UjX ■ VxUk + UkX ■ VxUj) + 0(e) 

jk 

(4.24) 



Note that 



□ 



X ' Vx(Ci • Vx^i) — ^ ^ Xj ^^{Ci^k'^XjXk 

j k 

Ci ■ Vx{x ■ VxUi) = Ci ■ VxUi + ^{Ci)k ^ XjUx^xk- 

k j 

Hence, the limit of the difference of the formulas (4.23), (4.24) yields 
{[F, G\[u])i = -2{F[u])i + {F[u])i = -{F[u])i. 

Corollary 4.1. Ifu{t,x) = {ui{t,x), . . . ,ui{t,x)) is a solution, then 

u^{t, x) = e^u{e^t, e^x) 

is another solution to the system of PDE (4-22). Hence, by performing the 
change of variables 

r = In 

^ = e-a=l!^ (4-25) 
. Wi = e~'^Ui = (1 - t)ui 
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blow-up to (4-22) occurs if a solution to 



(Wi)r = -(Ci + V) ■ ^riWi + ^ aijkWjWk - Wi 

jk 

approach a steady state as r goes to infinity. 
Proof. By Theorem 4.5 the vector field 

■ .9 „ d 

is the generator of symmetry associated to G, hence fixed 9, the graph of the 
solution is mapped by v in the graph of another solution: 

= {exp(^v)(f, x,u) : u = u{t, x)} 

that is 

t{e) = e'H 
x{e) = e~^x 
u{t{e),x{e)) = e^u{t,x) 

thus, the new solution is u^{t, x) = e^u{e^t, e^x). 

As far as the second part of the corollary is concerned, we shall use Theorem 
4.3. Suppose that there exists a solution w = {wi{t,t]), . . . ,wi{t,t])) to the 
system 

drW = F[w] - G[w] 

which approach a steady state as r ^ oo. Then the transformation by exp(Tv) 
with T = In is the graph of the solution 

u{t, x) = exp(TG)w(r, ri) 

which corresponds to the change of variables 

{X = e~'^r] = (1 — t)ri 
r Wi 
Ui = e'^Wi = J— ^ 

which yields (4.25). □ 



Chapter 5 



The two dimensional 
Broadwell model 

5.1 The discrete Boltzmann equation 

Consider the simplified model of a gas whose particles can have only finitely 
many speeds, say Ci,...,CAr e M". Call Ui = Ui{t,x) the density of particles 
with speed Cj. The evolution of these densities can then be described by a 
semilinear system of the form 

dtUi+Ci-'Vui = '^aijkUjUk i = l,...,N. (5.1) 

Here the coefficient aijk measures the rate at which new -i-particles are created, 
as a result of collisions between j- and /e-particles. In a realistic model, these 
coefficients must satisfy a set of identities, accounting for the conservation of 
mass, momentum and energy. 

Given a continuous, bounded initial data 

u,{0,x)=u,{x), (5.2) 

on a small time interval t G [0, T] a solution of the Cauchy problem can be con- 
structed by the method of characteristics. Indeed, since the system is semilinear, 
this solution is obtained as the fixed point of the integral transformation 

Ui{t,x) = Ui{x — Cit) + / ajjft MjUfc (s, X — Ci{t — s)) ds . 

For sufficiently small time intervals, the existence of a unique fixed point fol- 
lows from the contraction mapping principle, without any assumption on the 

constants aijk. 

If the initial data is suitably small, the solution remains uniformly bounded 
for all times [4]. For large initial data, on the other hand, the global existence 
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and stability of solutions is known only in the onc-dimcnsional case [3, 32, 45]. 
Since the right hand side has quadratic growth, it might happen that the solution 
blows up in finite time. Examples where the L°° norm of the solution becomes 
arbitrarily large as t ^ oo arc easy to construct [35] . In the present chapter we 
focus on the two-dimensional Broadwell model and examine the possibility that 
blow-up actually occurs in finite time. 

Since the equations (5.1) admit a natural symmetry group, one can perform 
an asymptotic rescaling of variables and ask whether there is a blow-up solution 
which, in the rescaled variables, converges to a steady state. This technique 
has been widely used to study blow-up singularities of reaction-diffusion equa- 
tions with super linear forcing terms [30, 31]. See also [36] for an example of 
self-similar blow-up for hyperbolic conservation laws. Our results show, how- 
ever, that for the two-dimensional Broadwell model no such self-similar blow-up 
solution exists. 

If blow-up occurs at a time T, our results imply that for times t — > T— one 

has 

II (,\\\ ^ 1 ln|ln(T-t)| 

hmL^ > 5 T-t ■ 

This means that the blow-up rate must be different from the natural growth 
rate = 0{1) ■ {T — t)~^ which would be obtained in case of a quadratic 

equation u = Cv?. 

In the final section of this chapter wc discuss a possible scenario for blow- 
up. The analysis highlights how carefully chosen should be the initial data, if 
blow-up is ever to happen. This suggests that finite time blow-up is a highly non- 
generic phenomenon, something one would not expect to encounter in numerical 
simulations. 

5.2 Coordinate rescaling 

In the following, we say that P* = {t* ,x*) is a blow-up point if 

limsup Ui{t,x) = oo 

for some i £ {1, . . . , N}. Define the constant 

C = max |ci| . 

i 

We say that {t*,x*) is a primary blow-up point if it is a blow-up point and the 
backward cone 

T={{t,x); \x-x*\<2C{t* -t)} 
does not contain any other blow-up point. 

Lemma 5.1. Let u ~ u{t,x) be a solution of the Cauchy problem (5.1)-(5.2) 
with continuous initial data. If no primary blow-up point exist, then u is con- 
tinuous on the whole domain [0,oo[xM". 
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Proof. If u is not continuous, it must be unbounded in the neighborhood of 
some point. Hence some blow-up point exists. Call B the set of such blow-up 
points. Define the function 

ip(x)= inf \t + C\x-^\} . 

By Ekeland's variational principle (see [1], p. 254), there exists a point x* such 
that 

(J 

V{x) >(p{x*)- — \x-x*\ 

for all X e M^. Then P* = ((^(a;*), a;*) is a primary blow-up point. □ 

Let now {t*,x*) be a primary blow-up point. One way to study the lo- 
cal asymptotic behavior of u is to rewrite the system in terms of the rescaled 
variables Wi = Wi{T, rj), defined by 

T = -\n{t*-t), 

rj = e-'x = (5.3) 
' _ t* -t 

Wi = e~'^Ui = {t* — t)ui. 

The corresponding system of evolution equations is 

drWi + {ci+r])-'VnWi= - Wi + '^aijkWjWk i = l,...,n. (5.4) 

Any nontrivial stationary or periodic solution w of (5.4) would yield a solution 
u of (5.1) which blows up at {t*,x*). On the other hand, the non-existence 
of such solutions for (5.4) would suggest that finite time blow-up for (5.1) is 
unlikely. 

5.3 The two-dimensional Broadwell model 

Consider a system on consisting of 4 types particles (fig. 5.1), with speeds 

Ci = (l,l), C2 = (l,-1), C3 = (-l,-l), C4 = (-l,l). 



The evolution equations are 

dtUi + C\ ■ Vui = Zt2U4 — M1M3 , 
^ dtU-i -h C3 • VU3 = - , /g 

dtU2 + C2 • V'U2 = Wl^S - W2U4 , 
dtU\ -\- C4 • VU4 = U\Uz — U2Ui . 

After renaming variables, the corresponding rescaled system (5.4) takes the form 

dtWi -I- (a; -I- l)dxWi + {y + l)dyWi = W2Wi — wiw^ — wi , 
^ dtWj, -I- (a; - l)dxW3 + {y - i)dyW3 = W2Wi - W1W3 -W3, 

dtW2 + {X + l)dxW2 + {y - l)dyW2 = WiWs - W2W4 - W2, 

dtWi -I- (a; - 1)9x^4 + (y + l)dyWi = W1W3 - W2'W4 - . 
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C4 



C3 



Cl 



C2 



Figure 5.1: Moving particles with prescribed speeds 



Our first result rules out the possibility of asymptotically self-similar blow-up 
solutions. A sharper estimate will be proved later. 

Theorem 5.1. The system (5.6) admits no nontrivial positive bounded solution 
which is constant or periodic in time. 



Proof. Assume 



< 'Wi{t,x,y) < K 



(5.7) 



for all t, x,y, i = 1, 2. Choose e = e ^'*/2, so that 



e < 



Define 



Qi4{t,y) = J ^[{ 



1 -ee2"^)wi(t,a;,y) + (l - ee-'^'^'')w4{t,x,y) 



dx . 



Qu{t) = sup Qu{t,y) , 
|y|<i 



Restricted to any horizontal moving line y = y{t) such that y = y+1 (fig. 5. 2), 
the equations (5.6) become 



dtWi + {x + l)dxWi = W2W4 - wiws - wi , 

dtW4 + {x — l)dxW4 = WiWs — W2W4 — W4 . 
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(-1,1) 




(1,1) 




(-1,-1) (1,-1) 
Figure 5.2: Interaction on a moving line 



A direct computation now yields 
di 



Qi4{t,y{t)) 



< 



e^'^^il + x)wi + e-^'^^(l - x)wi 



< —EK 



€^'^'^{1 + x)wi + e~^''^(l - X)W4 



- SK e ^"^(1 - x)'W4 dx 

pO rl ^ 

en / e^'^'^ {I + x)wi dx + I £Ke~'^'^'^Widx + I SKe'^^^'widx 
Jo J~i Jo 



< —EK 



e^''^(l + x)wi + e"^''^(l - x)wi 



dx . 



Call 



I{t,y) = j wx{t,x,y)dx. 
The definition of e and the bound (5.7) on w\ imply 



1 + x)wi dx > 2s J (1 + x)wi dx 

7-1+1 /k 
>2s {l + x)Kdx 
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Prom this, and a similar estimate for Wi, we obtain 



/: 



dx > — i widx] H — I / Widx 



2 / „l s 2 



1 \ ^ ^ / /•! 



K 2 



Since e < k ^, this yields 



^^Qu{t,y{t)) < ~Qli{t,y{t)) . (5.8) 



Observing that the Cauchy problem 

i = — -z"^, z{0) = 4:K 

has the solution 

by a comparison argument from (5.8) we deduce 

(3l4(t) <(^ + ^t 



Since 

»i ,.1 



4:K 2 



wi(t, a;, y) dxdy < AQult), 
and since a similar estimate can be performed for all components Wi, we conclude 

J J ^ Wi{t, X, y) dxdy < 4 (^^ + ■ (5.9) 

The right hand side of (5.9) approaches zero as t — !■ cx3. Therefore, nontrivial 
constant or time-periodic solutions of (5.6) cannot exist. □ 



5.4 Refined blow-up estimates 

If {t*,x*) is a blow-up point, our analysis has shown that in the rescalcd coor- 
dinates T, ^ the corresponding functions Wi must become unbounded as r — ^ oo. 
In this section we refine the previous result, establishing a lower bound for the 
rate at which such explosion takes place. 

Theorem 5.2. Let u be a continuous solution of the Broadwell system (5.3). 
Fix any point {t*,x*) and consider the corresponding rescaled variables r, ^, Wj. 
If 

max WiiT,£i,£2) < 9lnT z=l,2,3,4, 
l£iM€2|<i ^ '"^ ^ - 
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for some 9 <1/A and all r sufficiently large, then 

lim Wi(r,0 =0 i = 1,2,3,4, 

r— ^oo 

uniformly for ^ e in compact sets. Therefore {t* ,x*) is not a blow up point. 



Since Wi = {t* — t)ui and r = | \n{t* — t)\, the above implies 

Corollary 5.1. If{t*,x*) is a primary blow-up point, then 

t* — t 1 
limsup \u{t,x)\ ■ — -T > -. 
x^x*, t->t*- ln|ln(t* — fj| 4 

Proof of Theorem 5.2. 

Let Wi = Wi{t,x,y) provide a solution to the system (5.6), with 

0<Wi{t,x,y)<eint^k{t) (5.10) 

for all t > to and x,y £ [—1,1]. The proof will be given in two steps. First 
we show that the norm of the components Wi approaches zero as t ^ oo. 
Then we refine the estimates, and prove that also the L°° norm asymptotically 
vanishes. 

STEP 1: Integral estimates. Consider the function 



Qu{t,y) = J 



1 r/ e'^Ht)(^-i)\ / g-2fc(t)(a: + l) 

1 wi{t,x,y) + I 1 ) Wi{t,x,y) 



dx 



with k{t) as in (5.10). As in the proof of Theorem 5.1, let 1 1-^ y(f) be a solution 
to y = y + 1. Then 



d_ 
It 



Qu{t, y{t)) = j\-(x- l)fc'e2'=W("-i)u;i + {x + l)k'. 

j-l / g2/c(t)(a:-l)\ 

+ / ( 1 ^ j [ - (a; + l)wix + W2Wi - W1W3 - wi] dx 

Jl / g-2fe(t)(x+l)\ 

+ J (1 1 [ - {x - l)w4x + W1W3 - W2W4, - Wi] dx 



r-1 / „-2fe(t)(x+l 

-.0 — — 

To estimate the right hand side, we notice that 



rl / 2fe(t)(x-l)\ rl 

A= h \ [{1 + x)wix + wi]dx > k{t) (a; + l)e2'=(^-^)M;i rfa; 

fl / g-2fe(t)(x+l)\ ~U 

B = J M - j [ix - l)wix + Wi] dx > k{t) y (1 - a;)e-2'=(^+i)«;4 , 



. dx 

C = J {W1W3 — W2Wi) ( - — — ) dx 



< kit) 



g2fe(x-l) g-2fe(x+l) 

'2 2 

rl g2fc(x-l) ~ pO g-2fe(x+l)' 

/ Widx + k{t) / Wi dx . 

Jo 2 J_i 2 
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Therefore, 



^Qi4(t, y{t)) 



[-{x - l)A;'e2'=(*)(^-i)M;i + [x + l)A:'e-2'=(f)(^+i)^i;4] dx 
-A-B + C 

< k\t) f [(1 - x)e2'=(*)(--i)ti;i + (1 + :r)e-2'=(*)(-+i)u;4] dx 



(5.11) 



If /c(f) > 1/2, we claim that the following two inequalities hold: 

(1 _ ^ym{x-l) < ^ _ g2fc(t)(a=-l) ^ 

(1 + a;)e"^'=W(a:^+i) ^ _ g-2fe(t)(x+i) _ 
To prove the first inequality we need to show that 

hk{s) = 1 - e"'' + se^'=" > for ah s € [-2, 0] . 

This is clear because /ifc(O) — and 

h'k{s) = e^^'{l-2k + 2ks)<Q se[-2,0] 

if fc > 1/2. Hence hk{s) is positive for s € [—2,0], as claimed. The second 
inequality in (5.11) is proved similarly. 

When t>to = e^l^^^^ one has k{t) > \ and hence 



|Qi4(t, t/(i)) < fc'(i)Qi4-^ y J(l+a;)e2^W(--i)t«i+(l-x) 
Setting / = widx, we obtain 

/ (l + x)e2^W(^-i)M;idx > / (l + x)e-^'''^*^k(t)dx = e-^'^C'^TT- 

Using the above, and a similar estimate for the integral of W4, we obtain 
k{t) 



^ j [(1 + x)e2'=(*)(--i)«;i + (1 - 3;)e-2'=(*)(-+i)w,4]da; 



> 



> 



-4/c(t) 
g-4fe(t) 



1 N 2 



■u;4(ia; ) 



(5.12) 



Q?4. 



Calling 



(5i4(i) = max Qii{t,y) 

|3/|<1 



5.4 Refined blow-up estimates 
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from (5.12) we deduce 

^QuW < k'it)Qu{t) —Quitf. 

Recalling that k{t) = 9lnt for some < ^ < 1/4, the previous differential 
inequality can be written as 

j^Q..<jQ..-^Ql.. (5.13) 

Notice that Qi4{to,y{to)) < 2fc(fo), and define the constant 

Ao = max {2fc(to)4"''^ 8(1 - 30)} . 

Then the function 

z{t) = Aot""^-^ 

satisfies 

^z{t) >^z- z{to) > Quito, y{to)) ■ (5.14) 

Comparing (5.13) with (5.14) we conclude 

Quit) < zit) t>to. 
This implies the estimate 



i: 



Wiit, X, yo) dx < 2Qiiit) < 2Ao t 



AO -I 



for i > io, i G {1, 4} and any yo G !]• An entirely similar argument applied 
to Qi2-i Q23, ■ ■ ■ yields the estimates 

j Wiit,x,yo)dx<2Aot^^-\ J Wiit,xo,y)dy < 2Aot^^-K (5.15) 

for i = 1,2,3,4, XQ,yQ G [—1, 1] and t > to- 

STEP 2: Pointwisc estimates. Using the integral bounds (5.15), we now seek a 
uniform bound of the form 

Wiit,x,y)<Co (5.16) 
for some constant Co and all a;, y e [— 1, 1], f > 0. 

To prove (5.16), let 1 1— > xit), 1 1— > yit)) € [—1, 1] be solutions of 
x = x + l, y = y + 1. 

Call 

Ait) = / iwi +W4)it,x,yit)) dx . 

Jx{t) 
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Prom our previous estimates (5.15) it trivially follows 

A{t) < 4Aot^^-K 
The time derivative of A{t) is computed as 
dA 

— = -{xit) + l){wi+wi){t,x{t),y{t)) 



+ 



it) 



[dtwi + {y{t) + l)dyWi + dtW4 + {y{t) + l)dyW4\ dx 



= - {x{t) + 1) {wi + Wi) {t x{t), y{t)) 

- [wi + {x + l)dx'wi +W4 + {x- l)dxW4\ dx 

Jx{t) ^ 

= -{x{t) + l){wi + W4){t,x{t),y{t)) - {wi+W4)dx 

Jx{t) ^ 

-2wi{t,l,y{t)) + {x{t) + l)wi{t,x{t),y{t)) + / w^dx 

Jxit) 

-2w4{t,l,y{t)) + {x{t) - l)w4{t,x{t),y{t)) + Widx 

Jx(t) 

< [x{t)-l-{x{t) + l)\w4{t,x{t),y{t)) = -2w4{t,x{t),y{t)) . 
This implies 



W4{t,x{t),y{t)) < - 



IdA 
2~dt 



(5.17) 



The total derivative of wi along a characteristic line is now given by 

d / , , . .s. 1 f —dA 

— uii (t, x[t), y{t)) = W2W4 — W1W3 — wi < W2W4 — wi < -W2 I 

kify (-dA^ 



dt 



< -wi + 



Wl 



dt 



In turn, for t > to this yields the inequality 



wi{t,x{t),y{t)) <e 



-(t-ta) 



wi{to) + J e'^''>k{s){-A'{s))ds 
<e-^'-'°'>[wi{to)+A{to)Hto)] 
+e-(*-*o) [ A{s){e'-*°k{s)y ds. 



(5.18) 



The first term on the right hand side of (5.18) approaches zero exponentially 
fast. Concerning the second, we have 

e-(«-*o)^ A{s)e'-'°{k{s) + k'{s))ds< j\-^*-'^2Aos^^-^ (^e\ns+^^ ds . 

This also approaches zero as t — > 00. Repeating the same computations for all 
components, we conclude that for some time ti sufficiently large there holds 



X 1 
Wi{ti,x,y) < - 



for all x,y & [—1, 1] . 



(5.19) 



5.5 A tentative blow-up scenario 
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By continuity, the inequalities in (5.19) remain valid for all a;,y in a slightly 
larger square, say [— 1 — e, 1 + e]. For t > ti we now define 

M{t) = max ^^Wi{t,x,y)] i = l, 2,3,4, a;,y e [-1 - e*"*^e, 1 + e*"*^e]| . 

From the equations (5.6) and (5.19) it now follows 



-M{t)<-M{t) + M^{t)<^, M{t,)<-. 



M{t) < 



< e-"" • e*i for all T>ti. 



Ml 

Returning to the original variables Ui = e^Wi, this yields 

Ui < e*i 

in a whole neighborhood of the point P* = {t* ,x*). Hence P* is not a blow-up 
point. □ 

5.5 A tentative blow-up scenario 

For a solution of the rescaled equation (5.5), the total mass 

/•I /•! 4 



(*) = / / XI ^' y) dxdy 
J-iJ-U=i 



may well become unbounded as f — > oo. On the other hand, the one-dimensional 
integrals along horizontal or vertical segments decrease monotonically. Namely, 
if t y{t) satisfies y = y +1, then 

— J [wi{t, X, y{t)) + W4^{t, X, y{t))\ dx <0. 

Similarly, if i; = a; — 1, then 

J [mit, x{t),y) + W4{t, x{t),y)] dy<0. 

Analogous estimates hold for the sums W1+W2 and W2+W3. Therefore, a bound 
on the initial data 

Wi{0, X, y) e [0, M] for all x,y ^ [-1, 1] , 

yields uniform integral bounds on the line integrals of all components: 



j Wi{t,x,y) dx < iM , j 



Wi{t,x,y)dy < AM . 
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(-1,-1) (1,-1) 
Figure 5.3: A possible interaction between particles 



If finite time blow-up is to occur, the mass which is initially distributed along 
each horizontal or vertical segment must concentrate itself within a very small 
region, thus forming a narrow packet of particles with increasingly high density. 
A possible scenario is illustrated in fig. 5.3. A packet of 1-particles is initially 
located at Pi. In order to contribute to blow-up, this packet must remain 
within the unit square Q. At P2 these 1-particlcs interact with 3-particlcs and 
produce a packet of 4-particles. In turn, at P3 these interact with 2-particles and 
produce again a packet of 1-particles. After repeated interactions, the packet of 
alternatively 1- and 4-particles eventually enters within the smaller square Q' . 
After this time, it interacts with a packet of 2-particles at P5 (transforming it 
into a packet of 1-particles) and eventually exits from the domain Q. 

To help intuition, it is convenient to describe a packet as being "young" 
until it enters the smaller square Q', and "old" afterwards. To maintain a 
young packet inside Q, one needs the presence of old packets interacting with it 
near the points P2, P3, P4 . . . On the other hand, after it enters Q', our packet 
can in turn be used to hit another young packet, say at P5, and preventing it 
from leaving the domain Q. 

As t 00, the density of the packets must approach infinity. One thus 
expects that most of the mass will be concentrated along a finite number of one- 
dimensional curves. Say, the packet of alternatively 1- and 4-particles should be 
located along a moving curve 714 (t, 6), where ^ is a parameter along the curve. 
The time evolution of such a curve is of course governed by the equations 



5.5 A tentative blow-up scenario 
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d d 

-714 = CI or -714 = C4 

depending on whether 714 consists of 1- or 4-particles. The presence of 
interactions impose highly nonhnear constraints on these curves. For example, 
the interaction occurring in P5 at time t implies the crossing of the two curves 
714 and 712, namely 

7i4(i,^^) =7i2(t,^) = A 

for some parameter values 6,9. The complicated geometry of these curves re- 
sulting from the above constraints has not been analyzed. 
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